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We examine three trace formulas, the spectral shift function, Friedel sum 
rule and Levinson theorem in potential scattering theory. Although these 
formulas are closely related to each other, they have been often formulated 
in different settings so far. We propose a unified treatment for them. First 
we construct the spectral shift function in an alternative way, and then we 
prove that the spectral shift function thus constructed yields Friedel sum rule 
and Levinson theorem in arbitrary dimensions without assuming spherical 
symmetry of the potential. 
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1 Introduction 



In physical systems, universal nature often reflects the global, geometrical 
(topological) structure of the system. For example. Gauss's law in classical 

electromagnetism is a consequence of the geometrical structure of the three- 
dimensional Euclidean space. It states that the flux, $, of the electric field, 
E, through any closed surface, S, is proportional to the total charge, Q, 
enclosed by the surface S: 

$ = J E-da = Cg 

with a constant C. 

In this paper, we study an analogue to Gauss's law in scattering theory in 
quantum mechanics. Let us consider a metal with a single impurity at zero 
temperature. The impurity potential scatters the conduction electrons, and 
changes their charge distribution. For a fixed Fermi energy Ep, the "excess 
charge", Z{8,p), due to the impurity is defined to be the difference between 
the total numbers of levels in the Fermi sea with and without the impurity. 
Then the excess charge, Z(£f), equals the total phase shifts, 0(£f), of the 
scattering matrix, S{Ef), for the impurity potential: 

(1) 0{Ef) = ^ log detSi&F) = Z{Ef). 

This is known as Friedel sum rule [14] in solid state physics [25] . 

Since the excess charge, Z{8,f), is formally written in terms of the trace 
of the difference between the spectral projection operators with and without 
the impurity potential, it is closely related to the spectral shift function 
(SSF) which was initiated by Lifshitz [26]. We briefly describe the previous 
construction of the SSF. Let H and Hq be a pair of self-adjoint operators. 
Then the SSF, C(A), is defined as a function on M satisfying the following 
property: If / G C^(M), then 

Tr[f{H)-f{Ho)]= j nm\)dX. 

Here we note that this formula fixes ^(A) up to an additive constant. The 
SSF is known to exist if, for example, (H + i)"'^ — (Ho + i)""^ is a trace class 
operator with some m > (see, e.g., Birman-Yafaev [4] or Yafaev [34]). 
Formally, the SSF is written 

(2) a>^) = Tr[EH{X)-EHo{X)], 

where Ea{-) denotes the spectral projection of a self-adjoint operator A. 
(This formal expression (2) is nothing but the excess charge!) It is well- 
known, however, Eh{X) — EH^){X) is not necessarily in the trace class, even 
when the above assumption is satisfied. 
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As is well known, there are two standard constructions of the SSF (see 
also Pushnitski [31] and references therein for a more sophisticated repre- 
sentation of the SSF). The first one is due to Krein who defines the SSF as 
a locally function on M. This construction requires relatively weak as- 
sumptions, and the definition is global in A. However, the existence of ^(A) 
for a fixed A € M is not obvious in this construction. The other construction 
is to compute the difference of the spectral functions. Namely, under certain 
conditions, one can define 

for A in a "regular" energy region. This method is used widely in the semi- 
classical and microlocal study of the SSF (see Robert [33] and references 
therein). The advantage of this method is that one can study the behav- 
ior of .^(A) in detail locally in A. On the other hand, .^(A) is not defined 
globally in A, and the method requires slightly stronger assumption on the 
perturbation. 

We also remark that the behavior of finite- volume spectral shift functions 
for a large volume is studied in ref. [24, 17, 18, 15, 16]. In particular, under 
a certain condition, a sequence of finite-volume spectral shift functions is 
shown to converge to the SSF in the infinite- volume limit [17, 18, 15, 16]. 

1.1 Main Results 

We propose another construction of the spectral shift function, ^(A) = 
^{\;H,Ho), for a pair of Hamiltonians, H = —A -|- V and Hq = —A, 
on L^(R"). We assume that the potential V satisfies 

(3) \V{x)\ < C(x)-", xeR"" 

with some a > n + 3 and some C > 0, where we have written {x) := 
yl + |xp. The idea for our construction is to show the existence of the 
boundary value of the perturbation determinant directly using the stationary 
scattering theory. This is a variation of Krein's construction, but we can 
prove that ^(A) is defined for each A G (0, oo) \ app{H) and continuous in 
the same region. 

As an application, we consider Priedel sum rule. We first define the 
finite- volume excess charge, Zr(A), due to the impurity potential, V, by 

Zr{X) := Tr [^r{Eh{X) - Eh,{X))M > 

where i^nix) = '9{x/R) is a cutofi' function with a large R and with -d G 
Co°(M") satisfying = 1 in a neighborhood of x = 0. Then we can prove 

Z{X) := lim Zr{X) = ^(A) for A G (0, oo)\c7pp(iJ). 
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Namely, the excess charge, Z(X), in the infinite- volume limit is equal to the 
SSF, '^(A). On the other hand, the total phase shift, ^(A), for the scatter- 
ing matrix, S(\), is equal to ^(A) from the Birman-Krein formula. From 
these, we obtain that Priedel sum rule (1) holds for £i? G (0, oo)\(Tpp(i?) in 
arbitrary dimensions. 

We also show that the SSF thus obtained yields Levinson theorem: Let 
n > 3. Suppose that 

\V{x)\ < Cixy^, X G 

with some a > 2n and some C > 0. Then, we can prove 

^(0) = lim e(X) = N<+q, 
A-*-+0 

where N< is the number of the bound states of the Hamiltonian H with an 
eigenvalues less than or equal to zero; the number q is given as follows: In 
dimensions n = 3, 4, 




2, for n = 3; 
1, for n = 4 



if the zero energy is a resonance for H, and q = otherwise. In dimensions 
n > 5, q = because of the absence of a resonance [20]. 

The result in three dimensions recovers several variants of Levinson the- 
orem [30, 28, 12, 13, 9, 29]. Our normalization of the phase shift is not 
necessarily the same as the previous ones. As far as we know, the statement 
in the higher dimensions n > 4 is justified for the first time although it 
has been believed to be valid [6, 5]. We stress that our results of the total 
phase shift in one and two dimensions recover the previous results, too. We 
present their precise statements in Section 4 because they are more compli- 
cated than the results in dimensions n > 3. In consequence, the total phase 
shift in the zero energy limit is always quantized to an integer or a half-odd 
integer in arbitrary dimensions. Since the total phase shift equals the excess 
charge, this implies that the excess charge in the zero energy limit is quan- 
tized to the same value. We also stress that, since the previous approaches 
to Levinson theorem in dimensions n < 3 did not rely on the usefulness 
of the SSF, they suffered from the existence of exceptional points of the 
spectrum of the Hamiltonian H. In consequence, their statements are not 
necessarily expressed in a sophisticated form including all the exceptional 
cases. We refine the previous results in dimensions n < 3. 

The present paper is organized as follows: In Section 2, we first describe 
our method to construct the SSF in three or lower dimensions, and then 
extend it to higher dimensions. In Section 3, we prove that the SSF is equal 
to the excess charge. In Section 4, we compute the total phase shifts in the 
zero energy limit. 
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2 Construction of the Spectral Shift Function 

We construct the SSF for potential scattering theory. First, we describe our 
abstract scheme for the construction, and then prove the existence of the 
SSF. Consider a pair of Hamiltonians, 

H = Hq + V, Ho = -A onL2(M"). 

We suppose the potential V satisfies the bound (3) with a > n + 3. We 

may allow V to have some singularities, but assume that it is bounded for 
simplicity. By the invariance principle, we construct ^(A) as 

e(A; H, Ho) = -m + M)-'; {H + M)-\ {Hq + M)-') 

with some integer ^ > and a sufficiently large M > 0, where ^{\]A,Aq) 
denotes the SSF for a pair A and Aq. We recall the SSF is defined as 

^(A;^,^o) = - lim -Im logA^^ (z), 

where A^/^g(z) denotes the perturbation determinant defined by 

^A/A,{z) = det [{A - z)(Ao - z)-^] for z G C \ {<j{A) U ^7(^0))- 

It is easy to see Ayi/^„(2;) is well-defined if A — ^0 is trace class, and it is 
analytic in z. Moreover, 

^AMo(^)^AMo(^) = [{A - z)-' - {A^ - z)-^] , 

and C(A; A, .4o) = if A > sup(o-(^) U a{Ao)) [or if A < inf(c7(A) U (t{Aq))]. 
Hence we have an expression of the SSF: 

(4) ^(A; A, Aq) = lim -Ira 1 \(.A - w)'^ - (Aq - w)'^] dw, 

z->-\+iO TT J^^ 

where 7z denotes a contour in C+ := {z £ C \ Im 2; > 0} such that 72(0) = 
k > sup{a{A)Ua{Ao)) and 72(1) = z. Note that this expression is consistent 
with the formal formula (2) by virtue of Stone formula. 

2.1 Dimensions n < 3 

First, we prove the existence of ^(A) in dimensions n < 3. In the next 
section, we treat the case in dimensions n > 4. In this section, we set i = 1, 
namely, 

A = {H + M)-\ Ao = {Ho + M)-^ 

with a sufficiently large (fixed) Af > 0. Then it is well-known A — Aq € Ji, 
where 3p denotes the p-th trace ideal. Hence ^a/Ao{z) is well-defined, and 
the above definition applies. Now the key estimate of our construction is 
the following: We denote 

li{z) = (z + M)-^ = {z + M)-\ 
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Proposition 1. Let A G (0, oo) \ app{H). Then 

lim Tr [{A ~^i{z)r'-iAo- ^{z)r'] 

exists, and the limit is continuous in (0, oo) \app{H). 

Remark. We do not prove (or claim) {A— fi{X+iO))~^ — {Aq — iJ,{X+iO))~^ G 
Ji. We only prove the existence of the limit of the trace. 

Now combining Proposition 1 with the formula (4) , we have the following 
result on the SSF: 

Corollary 2. The SSF, ^(A), exists for A G (0,oo) \ app{H), and ^(A) is 
continuous in (0, oo) \app{H). 

Proof of Proposition 1. For z ^ (t{H) U a{HQ), 

{A - fi{z)r^ - {Ao - Kz))-^ = -{A - fi{z)r\A - Ao){Ao - Kz))'^ e Ji. 

Let 3/2 < P < {a — n)/2. Then by the standard commutator computations, 
we see 

(5) W := {x)'^{A - Ao){xf = -{x)^{H + M)-^V{Ho + M)-\xf G Ji. 
Then we write 

(6) Tr [{A -^^{z)r'-{Ao-Kz)r'] 

= -TV [{A - fi{z))-'{x)-^W{x)-^{Ao - fi{z))-' 

= -Tr [wixrf'iAo - fx{z))-\A - fi{z))-\x)-^' . 

Now in order to complete the proof, it suffices to show the following lemma. 

□ 

Lemma 3. For A G (0, oo) \ app{H), 

lun (xyf'iAo - n{z)r \A - ^i{z))-\x)-P 

exists in 5(L^(M")), and the limit is continuous in (0, oo) \ (jpp{H). 
Proof. We note 

^0 - = {Ho + M)-^ -{z + M)-i = -{z + M)-\Ho - z){Ho + M)"^ 
and hence 

{Ao - fi{z))-^ = -{z + M){Ho + M){Ho - z)''' 

= -{z + M)-{z + Mf{Ho - z)-^. 
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Similarly, we have 

{A - n{z))-^ = -{z + M)-{z + Mf{H - z)-^ 
= {z + M) - {z + Mf{Ho - z)-^ 

+ {z + M'f{Ho - z)-'^V{H - z)-^. 

Thus we have 

(^0 - ^l{z))-\A - n{z))-^ = ao{z) + ai{z){Ho - z)-^ + a2{z){Ho - z)-^ 
+ a3{z){Ho - z)-W{H - z)-^ + ai{z){Ho - z)-^V{H - z)-\ 

where aj{z) are polynomials in z. Since 

{xy\Ho- z)-^{x)-\ {x)-^{Ho- z)'^{x)-f^, and (x)'^ {H - z)-^x)-^ 

(with 7 > 1/2) are bounded and continuous in a complex neighborhood of 
A in C+ (see Agmon [1], Reed-Simon [32] Section XIII. 8), we conclude the 
assertion. □ 

2.2 Dimensions n> 4 

If n > 4, we set 

A={H + M)-\ Aq = {Ho + M)-^ 

with £ eZ such that n/2 -!<£<. n/2. Then we have 

e 

A-Ao = - + M)-^V{Ho + M)-^-^+^ 

i 

= - ^{Ho + M)-^V{Ho + M)-^-^+^ 

e. j 

+ + M)-^V{Ho + M)-^-i+V(ifo + M)- 

j=i k=i 



-e-i+j 



Iterating this procedure i!-times, and using the fact V{Hq + M)~^ G Jp for 
p > n/{2j), we learn A — Aq e Ji. Then the main part of the proof of 
Proposition 1 can be modified accordingly. 

In order to modify the proof of Lemma 3, we use 

e 

Ao - ii{z) = - Yiz + M)-^{Ho - z){Ho + M)-'-^+^ 
= -{Ho-z){Ho + M)-^Lo{z), 
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where ijl{z) = {z + M) ^, and we have written 

e 

Since Re {z + M) > M if z ~ A, Lq{z) is invertible, and we obtain 

= -L^\z) [l + {z + M){Ho-z)-^] . 

We also write 

i 

L{z) = Y^{z + M)-\H + M)-^+K 

Then we have 

{A^-ix{z)r\A-ix{z)r' 

= L^^{z)[ao{z) + ai{z){Ho - z)'^ + +a2{z){HQ - z)'^ 
+ a3{z){Ho - z)-^V{H - z)-^ + a4{z){Ho - z)-^V{H - z)-^jL-\z) 

with some polynomials aj{z) in z. Moreover, using the standard weight 
estimates, 

we can carry out the same argument as in the proof of Lemma 3. Conse- 
quently, we have: 

Proposition 4. Let A G (0, oo) \app{H). Then 

lim Tr[{A-ti{z)r'-{Ao-f,{z)r'] 

exists, and the limit is continuous in (0, oo) \app{H). Moreover, the SSF, 
^(A), exists for X G (0,oo)\app{H), and^{X) is continuous in{0,(x>)\app{H). 

3 Friedel Sum Rule 

In solid state physics [25], the difference of the number of the states given 
by the right-hand side of (2) has been often called the excess charge. In this 
section, we define the excess charge, and show that it is equivalent to the 
SSF. Besides, the SSF is equal to the total phase shift ^(A) which is given 
by 

e^'^^^W =det5(A), A>0, 
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where iS'(A) is the scattering matrix. By the invariance principle and the 
Birman-Krein formula, we have 

e{X) = a\;H,Ho) = -aK>^);A,Ao) 

with 6{X) = for A < cr{H). Therefore, the excess charge is equal to the 
total phase shift. This is nothing but Friedel sum rule. 

To begin with, we introduce a cutoff function "du^x) = 'd{x/R) with a 
large R > and with 'd G Co°(M") satisfying i? = 1 in a neighborhood of 
X = 0. Then the excess charge is defined by 

Z(A) := lim Tr [^r{Eh{X) - Eho{X))^r] , 
R-^oo 

where Ea{X) denotes the spectral projection: X(-oo,X]{^)- We want to show 
that the above limit exists, and that it is equivalent to the SSF under certain 
assumptions. 
We denote 

Zr{X) = Tr [MEh{X) - Eho{X))M 
for A > 0. Using the notation of Section 2, we recall 

EniX) = 1 - EAiniX)) = - lim Im- [ {A- w)-^dw. 
Hence we have 

Zr{X) = - lim Im - [ Tr [^r{{A - w)'^ - {Aq - wy^)^R]dw 
= Im- / Ti\^R{A-wy\x)-f^W{x)-^{Ao-w)-^^R]dw, 



since '&r{A — w) ^{x) ^ and (x) ^{Aq — w) ^'Or have norm limits as w 
/x(A) + zO. In particular, this implies the existence of Zr{X). We show 

Theorem 5. Let X G (0, oo) \ app{H). Then 

hm ZR{X) = aX;H,Ho). 

R-^oc 

Proof. From (4), (6) and the above representation of Zr{X), we have 

ZR{X)-aX;H,Ho) 
= ZR{X) + ^{fi{X);A,Ao) 

= lm- [ Ti[W{x)-^{Ao - wy\'di - 1){A - wy^xyl^jdw. 

Thus it suffices to show 

\\{xy^{Ao - wy^l - ^i)iA - wy^xy^W ^ O as i? ^ oo 
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uniformly in w E 7a- If w is away from a{A) U a{Ao), then 

||(x)-^(ylo - - ^i){A - w)-'{x)-^\\ 

= \\{x)-^iAo iv)-\x)^i{x)-^{l - i}i)){A - w)-\x)-^\\ 

< \\{x)-^{Ao - w)-\x)^\\ ■ |Kx)-^(l - -WiA- w)-m = 0{R-^) 

locally uniformly in w. Thus it suffice to consider the case w ~ /u(A) it iQ. 
Analogously to the argument in Section 2, it suffices to show 

\\{x)-f^{Ho- z)-^{l-'di){H - z)-^{x)-^\\ ^0 as iil^oo 

if z ~ A ± iO in C± = I ± Im z > 0}. Since 

{x)-P{Ho - z)-\l - .J^iH - z)-\x)-^ 
= {x)-^{Ho - z)-\l - i)i)iHo - z)-\x)-^ 

- {x)-^{Ho - z)-\l - ^i){Ho - z)-'V{H - z)-\x)-^ 

Theorem 5 now follows from the next lemma. □ 

Lemma 6. There exist U; a neighborhood of X±iQ in C±, e > and C > 
such that 

||(x)-^(i/o - z)-\l - ^i)iHo - zy^xr^W < CR-' 

for z G It. 

Proof. We consider the case A + zO only. The other case is similar. It is easy 
to observe that it suffices to show 

\\{xr^{Ho - zr\l - ^i)r,{Ho){Ho - < CR'^, 

where tj G C^((0,oo)) such that = 1 in a neighborhood of A. In order 
to show this, we use a Mourre-type microlocal resolvent estimate of Isozaki- 
Kitada [19]. See also [27]. We use the free case only, which is considerably 
simpler than their general results. Let p± G C°°([— 1, 1]) such that 

p+{t) + p-it) = 1; p±{t) = ii±t<~. 
We also set S and e so that 

We write 

= R^x)-^ (1 - wr{x)]'} p±{x ■ iuie), 
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where x = x/\x\. We quantize p± by the usual Kohn-Nirenberg pseudodif- 
ferential operator calculus: 

P±f{x)=p±{x,D,)f{x) = (27r)-"/2y p±{x,Oe^^<f{£,)d£,. 

Then we have 

(1 - di)rj{H^) = R-%xnP+ + P_), 

and p_(a;,^) satisfies the assumptions of Theorem 1 (or Theorem 1.2) of 
[19], uniformly m. R> 1, and we have 

\\{xYP-{Hq- z)-^{x)-^\\ <C forz~ A,zGC+. 

Similarly, we can apply the same argument to instead of P_. 

In fact, {{xY P-\-{x)~^)* is also a pseudodifferential operator, and its symbol 

can be computed by asymptotic expansions, up to an error of 0((a:)~°°). 
Let p+{x, ^) be the symbol of {{xY P-\-{x)~^)* . Then it has the same support 
property with p+(x,^). In particular, we can show 

P+{x,0 =P+{x,C) + i{9x ■ d^)p+{x,0 - id^p+{x,0 ■ {ex/ {xf) + 0{{x)~'^). 

Thus we obtain 

\\{x)-^{Ho- zy\xYP+{x)^~'\\ <C for z~A,zGC+, 

as well. Combining these, we have 

\\{x)-^{Ho - zr\l - ^i)r,{Ho){Ho - zr'{x)-^\\ 

= R-%x)-^{Ho - z)-'{xY{P+ + P-){Ho - z)-\x)-^ 

+ \\{x)-^{Ho - z)-\xYP+{x)^-'\\ ■ \\{x)-'-^-'HHo - 
< CR-' 

for z ~ A, z e C+. □ 

4 Levinson Theorem 

In this section, we compute the SSF in the zero energy limit, i.e., the total 
phase shift at the zero energy. As a result, we show that the total phase 
shift at the zero energy is quantized to an integer or a half-odd integer. Our 
results are summarized as the following three theorems: 
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Theorem 7. Let n = 1. Suppose that the potential V satisfies 



(7) < Const.(x)-", xGM", 

with some a > 5. Then, the total phase shift at the zero energy is given by 

0(0)= lim 0(A) = iV<-f, 

where is the number of the bound states of the Hamiltonian H with a 
strictly negative eigenvalue; q = if the zero energy is a resonance for H , 
and q = 1 otherwise. 

The proof is given in Appendix A. 

Remark . As is well known, there appears no bound state and at most a 
single resonance state at the zero energy in one dimension [7, 8]. The total 
phase shift 0(0) at the zero energy in one dimension was derived by [7, 8]. 
The present assumption (7) on the potential V is slightly different from 
theirs [8]. 

Theorem 8. Let n = 2. Suppose that the potential V satisfies 

(8) \V{x)\ < Const.(x)-", X G M", 
with some a > 6. Then, 

0(0) = lim 0(A) = N< + Nr, 

where N< is the number of the bound states of the Hamiltonian H with the 
eigenvalue less than or equal to zero, and is the number of the resonance 
states which show asymptotics, 

, C - X 

\x\^ 

for a large \x\. Here c is a two- component vector in . There appear at 
most two resonance states showing the above asymptotics because the number 
of linearly independent vectors is not beyond two in two dimensions. 

The proof is given in AppendixB. 

Remark . There appears another resonance state [6] which shows asymp- 
totics, V ~ Const., for a large |x|. Surprisingly, this state docs not con- 
tribute to the total phase shift at the zero energy [6]. A Levinson-type 
theorem in two dimensions was derived by [6]. In the case that the zero 
energy is neither an eigenvalue nor a resonance for the Hamiltonian H for a 
class of potentials, Cheney [10] proved Levinson theorem in two dimensions. 
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Theorem 9. Let n > 3. Suppose that 

(9) < Const.(x)-", xGM" 
with some a > 2n. Then, 

(10) ^(0) = lim ^(A) = N< + q, 

A-5>+0 

where N< is the number of the bound states of the Hamiltonian H with an 
eigenvalues less than or equal to zero; the number q is given as follows: In 
dimensions n = 3, 4, 

q=[h^ forn = 3; 
1 1, for n = 4 

if the zero energy is a resonance for H , and q = Q otherwise. In dimensions 
n>b,q = Q because of the absence of a resonance [20]. 

The proof is given in Sections 4.1, 4.2 and 4.3 in which we treat the cases 

with n = 3, n > 5 and n = 4, respectively. 

Remark . As is well known, there appears at most a single resonance state 
at the zero energy in three [28] or four dimensions [21]. Under an assump- 
tion on the existence of the trace of a time delay operator, a Levinson-type 
theorem was obtained by Osborn and Bolle [30], and BoUe and Wilk [9]. A 
relation between the phase shifts at the zero and the high energy limits was 
obtained by Newton [28, 29]. His derivation is based on the modified Pred- 
holm determinant. In the case that the zero energy is neither an eigenvalue 
nor a resonance for the Hamiltonian H for a class of potentials, Dreyfus 
[12, 13] proved that the total phase shift at the zero energy is equal to the 
number of the bound states with a strictly negative eigenvalue under a sim- 
ilar normalization condition to ours for the SSF. See also Theorem 1.14 in 
Chap. 9 of [35]. 

We also remark that, quite recently, Bellissard and Schulz-Baldes [3] 
obtain a Levinson-type theorem for tight-binding Hamiltonians with an im- 
purity on with n > 3. For related works, see the references in [3]. 

In order to compare Theorem 9 with the result obtained by Newton [28] , 
we need the high energy asymptotics of the total phase shift ^(A). Our 
expression of the SSF also yields the well known asymptotics: 

Theorem 10. Let n = 3. Suppose that in addition to the bound (9) in 
Theorem 9, the potential V satisfies 

(11) \d'^V{x)\ < Const. (x)-"' 

for all multi-indices k satisfying \k\ = 1,2 and some a' > 3. Then 

(12) ^(A) = 9{X) = / d'xV{x) + ^^-^ I d'xV{xf + o(A-V2) 
for a large A. 
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The proof of Theorem 10 is given in Appendix C. Since one has 



hm 

A— >-oo 







immediately from Theorem 10, the result (10) for n = 3 coincides with 
that obtained by Newton [28]. Thus the result (10) is nothing but Levinson 
theorem in three dimensions. The high energy asymptotics (12) should be 
compared to the result by Jensen [22]. Since he treated the derivative ^'(A) 
of the SSF, he did not show that the constant term in the expansion (12) 
is vanishing [35]. His assumptions on the potential V are slightly different 
from ours. For related works, see the references in [22] and [35]. 

4.1 Proof of Theorem 9 in the Case of n = 3 

As in the proof of Lemma 3, we have 

(^o-^)-^ = -(^ + M)+ ^" + y 

z — Ho 

and 

Substituting these into the expression (4) of SSF, we have 
(13) aX;H,Ho) 

= - lim - Im /" Tr \(z - H)'^ - (z - Hq)-^] (z + Mfdw 



lim -Im [ Tr \(z - H)'^ - (z - Hq)-^] 



dz, 



where we have used w = {z + M)~^. In order to compute the total phase 
shift ^(0) in the zero energy limit, we introduce the excess charge at the 
zero energy as 

(14) Z° = limlim— /" dzTr\(z- H)-^ - (z- Ho)-'^] , 

where the path 7^ ,5 is parametrized as 2; = ee^^ with S < <p <2Tr — S. Then 
the total phase shift 6{0) is given by 

^(0) = lime (A) = Z^ + N<, 

A4-0 

where is the number of the bound states of H with a strictly negative 
eigenvalue, A < 0. 
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The difference of the resolvents in the trace is written 

1111 

V- 



z-H z-Hn z-H z-H, 



^ .{z-Ho)-^V- ^ 



z — Hq z — H z — Hq 

Write 

(15) u = \V\^/'^s^V and v = \V\^/'^. 
Then, one has 

(16) {z - Ho){z - H)-^V = [1 + V{z - H)-'^] V 

= [l + uv{z - i?)"^] uv 
= u[l + v{z- H)~'^u] V 

for z ^ o-{H) U (t{Hq). These observations yield 

Tt[{z-H)-^-{z-H,)-^] 

= Tr {(z - Hq)-^u [1 + v{z - H)-^u] v{z - Hq)-^] 
= Tr {v{z - Ho)-^u [1 + v{z - H)-^u\ } . 

Here we have used the fact that the operators, {z — Hq)^^u and v{z 
are Hilbert-Schmidt class for z ^ a{HQ). Further, using 

(17) [1 + v{z - H)-^u] = [1 - v{z - Ho)-^u]-^ for z ^ a{Ho) U 
one has 

(18) Tr[{z-H)-'-{z-Ho)-'] 

= TV {v{z - Ho)~^u[l - v{z - Ho)~^u]-^} . 

Write ^ 

V —u = vRqU + vRqU, 

z — Ho 

where is the operator whose integral kernel is given by 

1 



Att \x\ 

with the characteristic function 



-Xr{x) 



1, if |x| < R; 
0, otherwise, 
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and the integral kernel of Rq is given by 



1 Q^V^l^l 
Ait \x\ 



Xr{x) 



i<f> 



with = 1 — Xr- We take R = e with 2/5 < < 1/2, recalling z = ee 
with a small e > and 5 < (I) <2it — 5. 

Lemma 11. We have vRqU G 32, and the Hilbert- Schmidt norm is bounded 



as 



\vR^u\\^ < Const. X £^+^' 



for a small £ with a positive constant S' . 
Proof. Note that 



<^\V{x)\xl{x-y)\V{y)\ 

= ^2^xf\V{x)\{yf\V{y)\{x)-\>{x-y){y)-^ 

Thus the statement follows from the assumption (9) for the potential V and 
R = e-i^ with fjL > 2/5. □ 

The integral kernel of the operator vRqU can be expanded as 



(19) 



4;^^(^) ix-y\ ^^^^ " y^^'^y^ 



47r ' \x — y\ ■' — ' ml 

' m=0 



1 °° 1 



47r \x — y 



Xr{x - y)u{y) - ^-^v{x)xr{x - y)u{y) 



+ —v{x)\x - y\xR{x - y)u{y) 



m=3 

Consider first the sum in the last line. Note that 
v{x){i^/z\x - y|)'""^X^(a; - y)u{y) 

= iV^{xf/Mx){yf/My){xr^^'\x - y\{yr'^^{iV^\x - y\r-\^{x - y) 
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for m > 3. The Hilbert-Schmidt norm of the operator of this right-hand 
side can be evaluated as 0(£^/^+'^ ) with some positive S', using 

{xr'\x + y\{yr' < {xr\\x\ + \y\){yr' < 2 

and e^/^i? = e^/'^~^ with fx < 1/2. Therefore the Hilbert-Schmidt norm of 
the operator corresponding to the sum in the right-hand side of the second 
equality of (19) is estimated as 0{e^^^ ). The first term in the right-hand 
side of the second equality of (19) is written 

- yHy) = 

In this right-hand side, the Hilbert-Schmidt norm of the operator corre- 
sponding to the second term can be estimated as 0{e^~^^ ) in the same way. 
Further the second and third terms in the right-hand side of the second 
equality of (19) are written 

- ^-^v{x)xr{x - y)u{y) + -^v{x)\x - y\xR{x - y)u{y) 
= -'^-^v{x)u{y) + -^v{x)\x - y\u{y) + 0{£^^^') 

in the sense of Hilbert-Schmidt norm. Combining these observations and 
the preceding lemma, we have 

(20) v—^u = -—vRI^\ - '-^vR^i\ + —vR^;^\ + 0(ei+^'), 

z — Hq Att Att Stt 

where are the operators whose integral kernel is given by Ixl"*"-*^ for 

m = 0,l,2,.... 

Let us consider the case when there are a single resonance state with 
the zero energy and Nq bound states with the zero energy because the other 
cases can be handled in much easier way. 

Let (fQj be the eigenvector of the operator uRo{0)v with the eigenvalue 
1 for J = 0, 1, ... , A^O) where Ro{z) = {z — Ho)~^. Namely, 

vRo{0)wpoj = ipoj, j = 0, 1, . . . , iVo. 

We take <^o,o to be the vector corresponding to the resonance state such that 

(u,¥?o,o) = J d^x u{x)(Po^q{x) ^ 0. 

Therefore the rest of ipoj for j = 1, . . . , TVq can be taken to satisfy 

(21) {u,ipo,j) = for j = l,2,...,iVo. 
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The left eigenvectors are given by 

<^o,i = (sgn V)(pQ^j for j = 0, 1, 2, . . . , Nq. 

As in the proof of Lemma 3.1 of [2], we can choose the vectors, <^o,j) for 
j = 0, 1, . . . , A^O) so that they satisfy 

{00,1,^0,3) = Sij. 

We define two projections as 

Po '■= ^Ofl{0o,o, ■■■) 

and 

No 

Pi :=^(^0j(^0j,---)- 

Further we write P = Pq + Pi and Q = 1 — P. Using (20) and the properties 
of the eigenvectors ^poj, we have 



(22) 



1 - V ^-rrU = f 1 - V ^"71"^ I <3 + ( 1 - V -T^^ ) P 

z- Ho \ z- Ho J V z- Hq J 

= (l- v^—uQ] Q + '4^vE!i\p - ^vRfKP 
\ z — Ho I Att ott 



+ 0(e 



1 - v^—uq] Q + '-^vR'^^uPo - ^vR^qKPo 
z — Ho J Att Stt 

_ A^i?(2)„p^+0(£l+*'). 
OTT 

We want to express the inverse of this operator in a similar expansion. Let 
^j^k be the No x A^o matrix which are defined by 

No 

(23) ^ ^i/(V'o/, V'Cfc) = Sj,k 

e=i 

with 

ipoj := Ro{0)u(po,j for j = 1, . . . , Aq. 

The functions ■0oj are all in L^(M^) because of the condition (21). (See 
[28, 2].) These functions are nothing but the eigenvectors at the zero energy 
for the Hamiltonian H. We note that, as in the proof of Lemma 3.1 of [2], 
it can be easily shown that the operator vRo{0)uQ has no eigenvalue +1. 



19 



Therefore the operator 1 — vRo(0)uQ is invertible. Relying on this fact, we 
introduce the following operators: 

with 



{u, <^o,o) 

OTT [U, (po,Q) 

and 

^4°^ = ^ E ^o,i*i/(<^o,^, • • • )vR^o^uil - vRo{0)uQ)-\ 



And 
with 



^(1) Airi 1 , ^ , 

-^1 := rT2<yCo,o(¥'o,o, • • • ) 



and 



■ No ^ 

'■=-^Y^ ^o,j-^jA'Po,e, ■ ■ ■ )vR^o^u-———^^o,o{ipofi, ■■■). 



Further, 



^^^^ •= I ^oj'^jAw---)- 

Using these operators, an approximate inverse for (22) can be written 

B = 5(0) +5(1)+ ^(2) _ 

Lemma 12. The following relation holds: 

(24) (^ipo J, vR^Q^ipo^k) = -Stt (V'oj, V'o.fc) for j,k = 1,2,..., Nq. 
Proof. Since 
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the integral kernel of Ro{zy is given by 



Using this and {u, (fioj) = for j = 1,2, . . . , Nq, one has 



(0o,j,vR^o^u(po^k) = -87rlim f <^o,i,^^ 



Ro{zf 



(1) 



-8tt lim {ipoj,vRo{zfuipo^k} 

-Stt lini {Ro{z)*v(po,j,Roiz)u(po,k) 

-Stt (V'oj,'0o,fe) 



□ 



From (20) and (22), one has 



1 



z-Hn 



(i-y^—uQ)Q + 0{e'/') 
\ z - Ho J 



= Q- 



-Q 



(2) fO,0(U, 



^ No 

No 



Q 
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Similarly, 



z-H, 







V z — Hq J Att 



- — 12 ¥'o,o(<^o,o, • • • )vR^qKq 



\{u,m,oW 
1 ^" 1 



1 



No 



+ StT fO,J^jA'fO,e, ■ ■ ■ )vR^o''u X, 2 ^0,0(^0,0, • • • )vR^o\Po 

foM^hmln < ^ it, \ d(2) y^o,o(^,---) ^ 

^ No 

+ ^0 + ^ </'o,,*i,£(<^o,£, • • • )vR^o^uPo + Oie'/' 
Further, we obtain 

V z- Ho 

+ 5(2) (-|-^;42)u)Pi+0(£^') 

iVo -j^ No 

= — X] <^o,j*j,^(<^o,^ • • • )vR^q\q - ^ X] <^o,j*j,K<^o,€, • • • )vR^qKPo 
j,e=i j,e=i 

^ No 

-g^Yl ^o,j^j4'Po,e, ■ ■ ■ )vR^o^uPi + 0{e'') 

No ^ No 

= -^Yl 'Poj^jA'Po/, ■ ■ ■ )vRfKQ - ^o,j'^jA'Po,e, • • • )vRfKPo 

j,e=i j,e=i 

+ Pi + 0{e^'), 

where we have used (23) and (24) for getting the last equality. From these 
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observations, we have 

b(i- V — K^u] =1 + D 
V z-Ho J 

with an operator D having the norm, 

\\Dh = 0{e''), 
with some positive constant S'. Therefore we obtain 

(25) (l-v — i — u] ={l + D)-^B = B-{l + D)-^DB. 

V z- Hq J 

Substituting this into the right-hand side of (18), we have 

(26) T^[{z-H)-^ -{z-Ho)-^] 

= Tr {vRo{zfu[B - (1 + D)-^DB]} 

= Tr vRo{zfuB - Tr {vRo{zfu{l + D)-^DB} . 

Consider first the first term in the right-hand side of the second equahty. It 
can be written 

Tr vRoizfuB = Tr vRo{zfuB^^^ + Tr vRQ{zfuB^^^ + Tr vRo{zfuB^'^\ 
Further, 

4 

Tr vRo{zfuB^^^ =^Ti: vRQ{zfuBf\ 

i=l 

In the same way as in the proof of (20), wc have 

(27) vRo{zfu = ^-^vR^o^u - ^vRf^u + 0{e^'). 

Combining this with the expression of Bf^\ one can notice that their con- 
tributions to the excess charge of (14) is vanishing except for that for Bf'\ 
Prom the definition of B^^ , we have 

Tr vRo{zfuBf^ = Tr vRQ{zfu (l - v — \r^Q] ■ 

\ z — Hq J 



Note that 



1 - V — ^——uQ 
z- Ho 



1 



1 - I 1 - V — 

z- Ho 



1 



-1 
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Therefore, 

TV vRo{zfuB[^^ 

= Tr vRoizfu + Tr vRQ{zfuv{z - HqY^uQ [l - v{z - HqY^uQ] . 

The contribution from the second term in the right-hand side to the excess 
charge can be shown to be vanishing by using the expansions, (20) and (27). 
The contribution from the first term is vanishing, too, because one has 

IV vM^fu = TV -J—V^— = -^-^ / d^yV{y) 

Z — Hq Z — Hq STTly/z J 

as in (88) in Appendix C. 

Next consider vRq{z)'^uB^^\ From the condition (21) and the expansion 
(27), the nonvanishing contribution to the excess charge comes from 



Stt^/z " 2z |(u,(/?o,o)r u , , 2z 

The corresponding contribution to the excess charge is 1/2. 

In the same way, the nonvanishing contribution for vRo{z)'^uB^'^^ comes 
from 

-— Tr vR^Kb^^^ = -— ^ TV vd^\^oj^j,i{^o/, . . . ) = -ATq, 
^ j,£=i ^ 

where we have used (23) and (24). This gives A'o for the excess charge. 

Consider the second term in the right-hand side of the second equahty of 
(26). Using the condition (21) and WDW^ = 0{e^ ), we can show, in the same 
way, that the corresponding contribution to the excess charge is vanishing. 

Combining these observations with (14) and (26), we obtain the excess 
charge, 

at the zero energy. Consequently, the total phase shift in the zero energy 
limit is given by 

(28) 0(0) = ^ + iVo + iV<. 

4.2 Proof of Theorem 9 in the Case of n > 5 

In order to prove Theorem 9 in the case of ra > 5, we choose j5 to satisfy 
n/2 < (3 < {a — n)/2. Then W of (5) is a trace class operator. Using the 
polar decomposition for W , we write 

A-Ao = {x)-^W{x)-^ = uv 
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with 

u={x)-^\W\^'^U and v = \W\^'^ {x)'^ . 

Here the unitary operator U satisfies U* = U because A — Aqis self-adjoint. 
The SSF (4) is expressed as 

^(A; A, Aq) = Im- [ dwTrUA- w)'^ - {Aq - w)~'^] 

= -Iin- [ dw Ti [{A - w)-\A - Ao)iAo-w)-^] 

= -Im - / dwTlcKA- w)~'^uv{Ao - w;)"^] 

= -Im - / dwTi[v{Ao-wy\A-wy^u]. 

One can easily check 

1 - v{A - wy^u = [1 - v{w - Ao)~^u]~-^ 

for w ^ <t{Aq) U a {A). This relation yields 

{Ao - w){A - w)-'^u =[1-{A- Ao){A - w)-'^]u 
= [1 — uv{A — w)~^]u 
= u[l — v{A — w)~^u] 
= u[l - v{w - Aq)~'^u]~^ 

for w ^ (t{Aq) U (t{A). By using this and the invariance principle, the SSF 
can be written 

^(A; H, Ho) =lm- f dw Tr {v{Ao - w)-^u[l - v{w - ^o)"^"]"^ } • 

First we want to rewrite {w — ^o)~^- Note that 

e 



w - A 







1 



z + Mj \Ho + M^ 
1 ^ 



^z + M Ho + M J j^J-^[(^{z + M) Ho + M\' 
where Q := exp[27ri/£]. Since one has 

1 ^__ ^ iHo + M)-C^{z + M) 

C^iz + M) Ho + M~ Cl{z + M){Ho + M) ' 
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the inverse is given by 
1 1 



H -1 



X^{z + M) Ho + M_ 



CKz + M). 



Ct{z + M) 



Ho + M 



{Ho + M) - C^{z + M) 



1 + 



{Ho + M)-CHz + M)\ 
Prom these observations, we have the expression, 

z + M 



(29) 



{w-Ao)-^ = {z + M) 



1 + 



Hq-z 



e-1 



k=l 



1 + 



CHz + M) 



{Hq + M)-Q{z + M)\ 



This yields 

(30) {xrP{w-Ao)-Hx)'P 

= {z + M) \{x)-^^ -{z + M){x)-^{z - Ho)-\x)-^ 



k=l 



-/3 



{Ho + M)-C^{z + M) 
Prom the definitions of u and v, we have 

(31) v{w - Ao)-^u = \W\^/^{x)-^{w - Ao)-\x)-^\W\^/^U. 

We write Gq{x; z) for the integral kernel of the resolvent Ro{z) = {z — 
Hq)~^. We decompose Ro{z) into two parts as Ro{z) = Rq{z) + Rq{z). 
Here the operator Rq{z) has the integral kernel Gq{x; z)Xfj^ix), and Rq{z) 
has the integral kernel Go{x; z)x^{x). In order to compute the excess charge 

at zero energy, we take z = ee^^ with a small £ > and 5 < (f) <2tt — 5 
with a small ^ > 0, and R = e~i^ with 1/3 < < 1/2. 

Using the integral representation of Go{x; z), one has 

\Gq{x-z)x^{x)\ < Const.e^+'^i, x G R", 
with some positive 5\. Combining this with /3 > n/2, we obtain 



{x)-^R^{x) 



0(e 



Using the expansion with respect to ^\x — y|, one has 

1 



{x)-^Go{x - y; z)x%{x - yM'^ = C,{x)-^- 



X — 2 



+ C2{X) 



-P. 



\x - y\ 



n—A 



xkix-y){vY 
XR{x-y){yY 



+ Ki{x,y), 
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where Ci and C2 are some constants, and Ki{x, y) is the integral kernel of a 
bounded operator. One can show that the norm of the operator for Ki[x, y) 
IS 0{e^+^^) with some positive 62- From these observations, we have 

{x)-^Ro{z){x)-^ = {xr^Ro{0){x)-^ + {x)-^^{0){xr^ ■ z + 0{e^^^') 

with some positive . Combining this, (30) and (31), we obtain 

(32) v{w - AqY^u 

= v[w{0) - Ao]-\ - v[w{0) - Ao]-'^u^iO) ■ z + 0{e^+^') 

= v[w{0) - Ao]-^u + j^vMO) - Ao]-^u + 0{e'+' ), 

where we have written w{z) = (z + M)~^. 

Let ifQ^j be the eigenvector of w[tf;(0) — Ao]~^^t with the eigenvalue 1 as 

v[w{{)) - Ao]~^u(fo^j = (fioj, j = 1,2,..., Nq. 

The left eigenvectors are given by 

<fO,j = U ifo,j, i = 1, 2, . . . , Nq. 

We choose these eigenvectors so that they satisfy 

We define 

No 

and write Q = 1 — Pi. Combining (32) with these properties of the eigen- 
vectors, we have 

(33) 1 - v{w - Ao)~'^u 

= [1 - v{w - Aor^uQ]Q + [1 - v{w - Aoy^u]Pi 

£z 

= |1 _ ^,[^^(0) - Ao]-\Q}Q - .j^v[w{0) - Ao]-\Q 

P7 

v\w{Q) - Aq]-'^uPi + 0(6^+^ ). 



In order to express the inverse of this operator in a similar expansion, we 
introduce 

V'Oj = [w{0) - Ao\~^Uipo,j, j = 1,2,..., Nq. 
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One can check that iJjqj G L'^(W^), and iI^qj are the eigenvectors of A with 
the eigenvalue w{0). This imphes HtpQj = 0. Let 'ifj^e be the A^o x -^o niatrix 
to satisfy 

Xl*i/(V'o,^V'o,fe) = <5j,fe- 

We introduce 
(34) 

B = {1 - v[w{0) - Ao]-^uQ}-^ 

No 

- ipoj^j^ifo/, ■ ■ ■ )v[w{0) - Ao]-^u{l - v[w{0) - Ao]-\Q}-^ 

Then, from (33), we have 

B[l - v{w - Ao)~\] 

No 

Prom the definitions of the projection Pi and of the vectors ipoj, one has 

No 

{'Po,e, ■ ■ ■ )v[w{0) - Ao]-'^uPi = Y{ipo,e,v[w{0) - Ao]-'^u(po,k){vo,k, ■■■) 

k=l 

No 

k=l 

Therefore, we have 

B[l - v{w - Ao)-^u] = Q + Pi + D = l + D, 

where D is a bounded operator with the norm of 0{e^ ). Consequently, the 
desired inverse is given by 

[1 - v{w - Ao)-^u]-^ = (1 + D)-^B = B-{1 + D)-^DB. 

Prom this and the expression (34) of the operator B, one notices that the 
nonvanishing contribution to the excess charge at the zero energy comes 
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from the second sum in the right-hand side of (34). Actually, we have 

1 f dz ^° 

= / V ^ Tr{v[wiO) - ^o]-'«¥'oj*i,^(^o,£, • • • )} 

1 f dz ^° 



from the definitions of the matrix and of the vectors tpQj . Adding this 
and the number of the eigenvalues below the zero energy, the desired result 
of Theorem 9 is obtained. 

4.3 Proof of Theorem 9 in the Case of n = 4 

In the present case, we choose 

A = {H + M)-^, Ao = iHo + M)-^ 

with a large M > 0. We also choose /3 to satisfy 2 < /? < (a - 4)/2. Then 
W of (5) is a trace class operator. From the condition /3 > 2, we can write 
P = 2 + S with a small S > 0. We define u and v in the same way as in 
Section 4.2. 

Since ^ = 2, we have 



{w - Ao)-^ = {z + Mf 

= {z + Mf 

_{z + Mf 
z-Ho 

from (29). This yields 
(36) {x)-^{w-Ao)-Hx)-^ 
= {z + Mf{x)-^ 

- {z + Mf{x)-^ 



1 - 



z + M 



z-Ho 

z + M 



z + M 



z + 2M + Hq 

z + M 
Z + 2M + Hq 



Z + 2M + Hq 
- 1 



- 1 



- 1 



z + M 



Z + 2M + Hn 



1 



z-H, 



-{x) 



-0 



{xY 



Z + M 
Z + 2M + Hq 



{x)-P-l 
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In order to obtain the expansion of {w — ^o) with respect to a smah 
\z\, we want to obtain the expansion of the resolvent Ro{z) = (z — Hq)~^. 
We write Gq{x; z) for the integral kernel of Ro{z). We decompose Ro{z) 
into two parts as Ro{z) = Rq{z) + Rq{z). Here, the operator Rq{z) has 
the integral kernel G{){x; z)x^{x) , and Rq{z) has Gq{x] z)x^{x). We take 
z = ee^^ with a small £ > and 5 < (f) <2it — 5 with a small ^ > 0. We also 
take R = e^^ with 

3 1 
6 + 5 2 

Lemma 13. We have 

{x)~^R^{z){xy^\l^^ < Const.e^+''' 

with a small S' > 0. 

Proof. Prom the integral representation of the integral kernel Go{x;z), one 
has 

\Go{x; z)xr{x)\ < Const.£3'^-^/2_ 
Combining this with ^ = 2 + 5, we have 

{xr^\Go{x-y;z)x^{x-y)\{yr^ 



< Const.£^'*-^/2^a 



{y)~^-'I^R-'"\ 



This yields 



< Const.£3^-V2+W2. 



{x)-PR>{z){x)-P 

On the other hand, we have 

3/x- l/2 + ^At/2 > 1 

from the condition fi > 3/(6 + 5). Combining this with the above bound 
yields the desired result. □ 

From the expansion of Go{x; z) with respect to v'^l^l) we have 

1 1 



Go{x; z)xr{x) = 



^zlog^ + Acz + ^loglxl 



47r2 



+ 0{e 



1+5'. 



with the constant, 



1 

8^ 



TT , 1 

2.+7-log2-- 
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and with a small 5' . Combining this with the above Lemma, we obtain 



+ ^{^r^R^oH^r^ + {x)-^^-\x)-^ + 0{e'+''), 

where the operators, i?Q^^ and i?o^\ have the integral kernels, 1 and log \x\, 
respectively. Therefore, from this expansion, (31) and (36), we have 

(37) v{w - Aoy ^u = v{w{0) - Ao]~'^u + z\og zAo + zAq + 0{e^^^'), 

where 

1 — -zr-r^ 7^ U 



2M + Ho_ 

and yio is a bounded operator. This yields 

(38) 1-viw- Aq)-^u = 1 - v[w{0) - Aq]-\ + Ki + Di, 

where 

Ki = — zlogz^lo — zAq 

and the operator Di satisfies ||i^i||2 = 0{e^~^^'). 

In order to obtain the inverse of 1 — v{w — Ao)~^u, let us consider the 
eigenvalue problem for the operator v[w{0)—Ao]~^u. Let ip be an eigenvector 
with the eigenvalue 1, i.e., 

'y['u;(0) — Ao]~^u(p = (p. 

The corresponding left-eigenvector is given by <^ = Uip. We set ip = [w{0) — 
Ao]~^u(p. Then, 

Ail^ = {[Ao- w{0) + w(0)] + uv} [w{0) - Ao]~^u(p 
= —uip + w{0)[w{0) — Aol^^uip + Uip 
= wiO)ij, 

where we have used A = Aq + {A — Aq) = Aq + uv. This implies 

{H + M)-^ip = M-'^tp. 

From the expression of {w — ^o)~^! one notices that, if (1, uip) ^ 0, then 
■i/j <^ L^(M^). This state ■0 is nothing but the resonance state at the zero 
energy. 

On the other hand, when {I, uip) = 0, ijj e L^(M^). This V' is a bound 
state at the zero energy. In order to see this fact, we note the following: 
The integral kernel of is given by 

1 1 
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For this integral kernel, we have 
1 1 



\x — 1 + \x\ 



^ 2|a:|H + l + H^ 
~ \x — + 

2|x||y| 



(1 + |xP)|x — yp |x — yP 



where we have used the inequality, 

2\x\ 



< 



l + - l + |x|' 

for getting the last inequality. Prom Lemma 2.3 of [20], both of two oper- 
ators, {x)~^A~^{x)~^^'~'^'' and {x)~^ {x)~^^~^\ are abounded operator. 
Combining these observations with the condition (l,u(p) = 0, one can show 
that V = [w{0) - Ao]-^u(p G L^{R'^). 

Thus, there appears at most a single resonance state at the zero energy 
in four dimensions [21]. In the following, we treat the case that there appear 
a single resonance state and N bound states at the zero energy because the 
rest of the cases can be handled in the same way. 

Let (foj, j = 0,1, . . . ,N, satisfy the following conditions: 

v[w{0) - Ao]~^uipoj = (fioj, 
{l,uipo.o) ^ 0, 

and 

(39) {l,u^oj) = fori = l,...,iV. 

The corresponding left-eigenvectors are given by 

and 

^o,j = ^o,j U ifoj , j = 1,...,N, 
where 'Nqj are normalization constants which wc choose to satisfy 

{^o,i, <fo,j) = i,j = 0,1,... ,N. 



We define 
and 



Po = </'o,o(<^o,o, ■ 



N 

and define P = Pq + Pi and Q = 1 - P. 
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Lemma 14. The following relation is valid: 

(<^o,fe,^o</fo,^) = -j^ii'o,k,i^o,e) fork,e = l,...,N. 
Proof. Prom the expansion (37), one has 

^{'fo,k, {viw - Ao)~^u - v[w{0) - Ao]~-^u} (po,i) 

foi k,£ = 1, . . . , N, where we have used the conditions (39). The left-hand 
side is computed as 

^(<^o,fe, {v{w - Ao)~'^u - v[w{0) - Ao]~'^u} (po^e) 

w(Q) — w 1 < 1-1 

^ -['Po,k,v{w - Aq) [w{0) - Ao\ uipo.e) 



z 

~ ^" —1 

K,fe((^* - ^o) uipo,k, iko,e) 



z 

In the hmit z ^ 0, this right-hand side goes to 

-^(V'o,fe,^o/)- 



Relying on this lemma, one can find N x N matrix $ such that 

N 

(40) ^ ^jk{'Po,k,'^ofo,e) = Sji 

k=l 

foi j,£ = l,...,N. 

Now, as an approximate inverse of the operator (38), we introduce 

5 = (1 - UKi) (1 - v[w{0) - Ao]~^uQy^ + n 

with 



1 ^ 

^0 - - XI fO,j^jk{'PO,k, ■■■)■ 



, 32vr-^ ^ ^ 

zlogzM^y^l,\il,u^o,o)r'' z.^^^^ 

Using the conditions (39) and Lemma 14, we have 

UKiP = P + E + PoxOi\loge\-^) 
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□ 



with 



N 



j,k=l 

Combining this with the expression of the operator (38) , we obtain 
B[l-v{w - Ao)~^u] 

= (1 - UKi)Q + (1 - UKi) {1 - v[wiO) - Ao]-^uQy ^ Ki 

= Q + UKiP + 0{£^') 

= 1 + H + Po xO(|loge|-^) + 0(£^'). 



In consequence, 



where 



[1 - v{w - Aq) ^u] 



1„.1-1 



1 



1 + D2 



:i-E)B, 



D2 = P X 0{\loge\-^) + 0{e^'). 

Let us calculate the excess charge at the zero energy. In the same 
way as the expansion of (37), we have 

d r A \-i dw . . ,_2 

—v[w — ^0) u = — —viw — Aq) u 
dz dz 

= {l + logz)Ao+Ao + 0{e^'). 

From the expression (35) for the excess charge Z^, it is enough to calculate 

1 



/ 



- / dzTr 



(1 + log^;)Ao +^0 



I + D2 



{1-E)B. 



Further, from the expression of the operator S, it is enough to calculate the 
following two contributions: 



1 



327r^ 



zM3K,o|(l,ii^o,o) 



rTr^lr 



1 



I + D2 



(1 - S)Po 



and 



Tr 



{l + logz)Ao + Ao 



1 + A 



il-E) 



N 



j,k=l 



except for the integral about z. The former is calculated as 



- + 0(£-l) 
Z 
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in the same way as in the above. Therefore, this leads to the excess charge 1 
from the resonance state. Using the conditions (39) and the definition (40) 
of the matrix the contribution from the latter becomes 

- + o(s \ 

This implies that the A'' bound states at the zero energy lead to the excess 
charge iV. 

A Proof of Theorem 7 

Since a > 5, one can find a small 5 > such that a > 5 + (5. As in 
Section 4.1, we take z = ee**^ with S < (f) < 2Tr — S, and R = with 
2/(4 + 5) <fi< 1/2. 

The integral kernel of the resolvent {z — Hq)~^ is given by 

Go(x;z) = -i^e^^l-l. 
2^^/z 

We decompose the resolvent Ro{z) = {z — Ho)~^ as 

Ro{z)=R>iz) + R^{z), 

where the operators, Rq{z) and Rq{z), have the integral kernels, 

Go{x; z)x^{x) and Go(x; z)x^(x), respectively. We write u = \V\^^^sgnV 

and V = |F|V2. 

Lemma 15. We have vRq{z)u G J2, and the Hilbert- Schmidt norm is 
bounded as 

\\vR^{z)u\\^ < Const. X e^/^+i' 
for a small e with a positive constant 5' . 
Proof. We have 

\Vix)\'/^^e^-^'\--y\xU^ - y)\V{y)\'/hgnV{y) 
<^jV{x)\\V{y)\x^{x-y) 

< l|y(x)Kx)^+^(x)-^-V(2/)l(y)'+'(y)-'-'x^(x - y) 
1 / 2 \ 

Since a — 4 — ^ > 1 from the assumption, one has 

\\vR>{z)u\\^ < Const. X £-'/^R-^-^/^ = Const, x e-i/2+(2+5/2),._ 
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The assumption, /x > 2/(4 + yields 



-l/2 + (2 + V2)/x> ^. 

Therefore the desired bound holds. 

As to Rq{z), we expand the integral kernel as 

J_e'^^''-y^XR{x - y) 

1 °° 1 

V ™— Q 

1 1 

1 + zv^|x -y\ - -z\x - + — ^(^v^|a; - 
z m! 



□ 



Note that 



rra=3 



oo oo 
* m=3 ^ m=3 

oo ^ 



m=3 



and 



<\\V{x)\^l\x)\yf\V{ytl\ 
Prom these observations, the integral kernel of vRq{z)u is written 



1 1 \/z 

-^—^v{x)u{y) + - - —v{x)\x - yfu{y) 

+ 0{e 



xkix-y) 



1/2+5' ^ 



in the sense of Hilbert-Schmidt norm. Using the same argument as in the 
above, we can replace the cutoff x^(a; — y) with the identity. Combining 
this with Lemma 15, we obtain 

(41) vRo{z)u = -^vR^^'^u + \vRf'^u - ^vR^^u + 0{e^l'^+^'), 

where i?Q™^ are the operators whose integral kernel is given by Ixp"^. 
Consider first the case with 



y" dxV{x) ^ 0. 
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Prom (14) and (18), the excess charge at the zero energy is given by 
(42) Z° = hmUm ^ [ dzTr {v(z - Ho)~^u[l - viz - Ho)-^u]-H . 

£4,0 54,0 ZTTZ Jj^ g 



We want to obtain the expression of the inverse of 1 — v{z — Hq) in the 
right-hand side. 
Note that 

l-v{z- Hq)-^u = 1 + li^p - M, 



where we have written 



2v^ 



P = {u,v) ^v{u,---), 



and 

(43) 



2 u 4i " 



The inverse of the right-hand side is formahy written 

-1 



1 + ^^p-m" 



2v^ 



i{u, v)P 
2y/z + i{u,v) 







M 


ri 



1 - 



2y/z + i{u,v) 



Note that 



1 - 



i(ii, v)P 
2-v/z -I- i{u,v) 



= 1-P + 



2y/z + i{u,v) i{u,v) 
_2^ + i{u,v) 2^ + i{u,v) 



2^Jz + i[u,v) 
where Q = 1 — P. Therefore, we have 
(44) [1 - v{z - Ho)-^u]-^ 



1-QM- 



2^z 



Since we have 
\-QM 



2^/z + i{u,v) 
2^ 



PM 



Q + 



2^z 



2^ + i{u,v)' 



2yfz + i(tt, V 



-PM 



= (1 - QMP) 



1 - QMQ - (1 + QMP) 



2^ 



2-v/z + %{u,v) 



PM 
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we obtain 



(45) 



1 - QM 



PM 



+ i{u,v) 
1 - QMQ - (1 + QMP) 



Note that 
(46) 



1 - QMQ - (1 + QMP) 



PM 



(l + QMP). 



+ i{u,v) 



PM 



1 - QMoQ + :^QvRfuQ 



2Vi 

i{u, v) 

where we have written 



1 + 



QMoP)PMo + 0(£V2+5')^ 



Mn 



1 



In order to find the inverse of the operator of (46), consider the eigenvector 
with the eigenvalue +1 for the operator QMqQ. The following lemma is due 
to BoUe, Gesztesy and Wilk [8]: 

Lemma 16. // the operator QMqQ has the eigenvalue +1, then the eigen- 
value must be simple and H correspondingly shows a zero- energy resonance. 

Proof. Let ipo be the eigenvector of QMqQ with the eigenvalue +1, and let 

where is the inverse of the Laplacian A. The integral kernel of 

is given by \x\/2. Clearly, Mq = vA~^u. Prom QMoQifQ = ipo, we have 

P(fo = and {u, ipo) = 0. Using these, we have 

QMoQifo = (1 - P)Moipo 

= Mo(po - PMQipQ 

= Mo(po - {u,vy^{u,Mo(po)v = (pQ. 

Therefore, 

(47) vtp = -ipo. 

Substituting this into the expression of the vector ijj, one has 
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Immediately, A^; = Vip, i.e., Hip = 0. Further, one obtains {u,(po) = 
— {V,ip) = and {xu,ipo) = —{xVjip) from (47). Using these relations, we 
have 

A- V^(x') = - / dy{y - x')V{y)iP{y) + - / dy{x' - y)V{y)iP{y) 

^ Jx' ^ J-oo 

= -x dyV{y)ip{y)+ I dyyV{y)tlj{y) + -{xu,(po) 

Jx' Jx' ^ 



X' J X 

OO -j^ 



dy{y - x')V{y)ip{y) + -{xu, ipo). 



Similarly, 

A-W'(x') = dy{x' -y)V{y)1P{y)-^{xu,^o)■ 



These yield ipeL'^i 

In order to show that ipo is unique, let ipg be another eigenvector with 
the eigenvalue +1, and let ip' be the corresponding resonance state of H. 
We choose their normalization so that they satisfy 

{u,v)~'^{u,Moi<fo - v'o)) + ^{xu, {(fo - (fo)) = 0. 
Then, we have 

POO 

iV(x') - < / dy\y - - ip'iy)\ 

Jx' 
roo 

< / dyi\y\+R)\V{y)my)-i;'iy)\ 

Jx' 

loi —R < x' < OO with a large R > 0. Iterating this, we obtain 

poo poo 

|^(x')-V'V)l< / dy^i\yi\+R)\Viyi)\ dy2{\y2\ + R)\V iy2)\ 

J x' Jyi 

dym{\ym\ + R)\V{ym)\\ip{ym) - ip'{ym)\ 

'Vm-l 



I 

Jy, 



< 



ml 



dy{\y\+R)\V{y)\ 



This implies tp = tp' . Hence, ipo = ip'^. □ 

In the following, we will treat only the case that H shows a zero-energy 
resonance because the other case can be handled in a much easier way. 
Let ipo be the eigenvector of QMqQ with the eigenvalue +1, i.e., 

QMoQ(fo = (fo- 
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Then the left eigenvector is given by (f>Q = 3\fo(sgny)(/?o- We choose the 
normahzation J^q so that ((^O)V'o) = 1- Wc define the projection Pq = 
'^o('^O) ■ ■ ■ )' Qo = 1 ~ -fb- Clearly we have 



(u, (/?o) = 



and 



PPo = PoP = 0. 
Prom these observations, we can write 



(48) 



where 



1-QMQ-{1 + QMP)- 



+ %{u^ V 

(1 - QMoQQo)Qo + V^Ki + Di, 



-PM 



2i 
{u,v)' 



Ki = -QvRfuQ + j^PMq + j^QMqPMq, 



2i 
{u,v) 



and the operator Di satisfies 

Pl||2 = 0(£V2+5')_ 

We introduce the approximate inverse, 



(49) 



P = (1 + K2){1 - QMoQQo)'^ + i^Po, 



for the operator of (48), where 



(50) 
and 

Note that 



N= (K^ 



{u, v) ^ I (m, Mo¥?o) 1^ + 7 1 {xu, (/?o) P 



1 -1 



K2 = ^PqvR^^^Q 



{u,v) 



PoMoPMo. 



2i 



PoKi = Po - K2Q0. 



Combining this, (48) and (49), we have 



(51) B 

where 

Do = 



1 - QMQ - (1 + QMP) 



2^ 



2^/z + i{u,v) 



PM 



I + D2, 



2U/z 



PqDi + ^{1 + K2)(l - QMoQQo)-'Ki + C»(£i/2+5'). 
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Therefore we obtain 



l-QMQ-il + QMP)- 



Combining this, (44) and (45), we have 

1 



PM 



-i„,i-i 



[1 - v{z - Ho)-'u] 



-B{1 + QMP) 



Q + 



I + D2 

2^ 



-B. 



P 



2^/z + i{u,v) 



I + D2' 

Let us compute the excess charge of (42) at the zero energy. Prom 

the above result, we have 

(52) TV {v{z - Hq)-'^u[1 - v{z - Ho)-^u]-^] 



1 



-BQ 



+ 



2^/z + ^{u,v) 



1 



I + D2 



B{P + QMP) \ . 



The first term in the right-hand side is written 



(53) 



Tvlviz-H, 



1 



oj u 



I + D2 



BQ 



= Ti:{v{z- Ho)~'^u 



1 



I + D2 



{1 + K2){1 - QMoQQo)-^Q 



+ 



Ti\v{z- Ho)~^u 



1 



I + D2 



2i^/z' 

by using (49) and PqP = 0. 

First, we will show that the first term in the right-hand side of (53) yields 
the vanishing contribution to the excess charge. To begin with, we note that 

we have 

(54) v{z - Ho)-^u = -'-z-'/\u, v)P - lz-^'\Rfu + 0{e-^/^+'') 

4 8 

in the same way as in (41). Formally substituting this expansion into the 
term, we obtain the desired result. But all of the terms so obtained cannot 
be expressed into a trace of a product of two Hilbert-Schmidt operators. By 
using the expansion and 

1 1 

I -Do 



I + D2 ' n + £>2 ' 
the first term in the right-hand side of (53) is written 



Tr <{ u(z - HoY'^u \l-D; 



1 



+ K2-D2 



I + D2 I + D2 

(1 - QMoQQo)-^q] 

Tr ^v{z - Ho)-M^ - QMoQQor^Q] + 0(£-V2). 
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Further, using 

(1 - QMqQQo)-^ = 1 + (1 - QMqQQqY^QMqQQq, 

we obtain 

Tr i^v{z - Ho)-'^u{l - QMoQQo)-^q] 

= Tr \^v{z - Hoy^uQ^ + 0{e^^/^) 

= Tr v{z - Ho)-^u - Tr v{z - HoY^uP + 0{e-^/^) 

= -'-z-'/\u,v) + ^z-'/\u,v)TrP + 0{e-^/^) = 0(£-V2). 

Thus the contribution is vanishing. 

Similarly, the second term in the right-hand side of (53) is computed as 

(55) -^Tt\v{z- Ho)-\- ^ 



= ^{(po,xv){xu,ipo) + 0(e"-^+*') 

where we have used PqP = 0, {(po, v) = {u, ipo) = 

Next, we treat the second term in the right-hand of (52). In the same 
way, we have 



1 / AK'')t^,p. 



X 



2^/z + ^{u,v)\ 2 J z \ I + D2 
(1 + K2){1 - QMoQQo)-^ + TT-^Po 
1 

2z %/z + 



(P + QMP)| +0(£-V2) 



-i{u,v) 4^3/2 \ 1 + ^2 J ^ 



-1+5' 



where we have used PPq = 0, PK2 = and P(l - QMqQQq)-^ = P. The 
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trace of the operator in the second term in the last hne is computed as 



Tr P- 



I + D2 



PqMP = Tr P 



I -Do 



I + D2 



PqM 



-Tr PD2- 



1 



-PqM 



I + D2 

= -v^Tr PK^-1—PqM + 0(£V2+5') 
i + U2 

= -2i{u, v)-^^/^Tr PMqPqMq + 0{e^l^+^') 
Substituting this into (56) , the value becomes 

Combining this, (50) and (55), all of the contributions turn out to cancel 
each other out. Thus the excess charge at the zero energy is vanishing. 
Next consider the case with 



/ 



dxV{x) = 0. 



We write 



P = v{u, ■■■) 
with u = |y|^/^sgny and v = Since we have 



v{z — Hq) 



P + M, 



we obtain 



l-v{z- Hq)~^u = 1 + -^P - M 

z 



= 1 + 



P 



1 - 



2v^ 



P M 



Here M is given by (43). Therefore, 



[1 - v{z - Ho)-\]-^ 



1 -1 



M + 



2v^ 



PM 



2^ 



P 



Since 



{u, Mv) = {u, vA-^uv) + 0(£^/2) = {V, A-^V) + 0{e 



l/2^ 



the matrix element {u, Mv) is nonvanishing for a small e. Relying on this, 
we define 

P = {u,Mv)-^PM, 
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and 



Q = l-P. 



Then, one can easily show that = P and QP = PQ = 0. From these 
observations, we have 



1 -M + 



PM 



1 



2yi 
i{u,Mv)P 







M 





1 



2^/z + i{u,Mv) 
in the same way as in the preceding case. Note that 

i{u,Mv)P ~ 2v^ 



i{u,Mv)P 
2^/z + i{u, Mv) 



1 - 



2v^ + i{u, Mv) 



Q + 



These yield 



1- 
Q + 



Q + 



2^/z + i{u,Mv) 

2^z 



P. 



2^ + i{u, Mv) 
2^ 



P 



M 



1 - QM 



2y^ + i{u,Mv) 
2^z 



P 



P 



2^ + i(u, Mv) 



-P + 



2^z 
-\ -1 

PM 

2v^ 



2^/z + i{u,Mv) 2^/z + i{u,Mv) 
where we have used QP = and PP = P. Note that 

2Vi 



P 



1 - QM 



2^/z + ^{u,Mv) 



PM 



(1 - QMP) 



1 - QMQ - (1 + QMP) 



2^ 



2^ + i(n, Mv) 



PM 



Therefore, the inverse is 



1-QM 



2^z 



2^/z + i{u,Mv) 
1 - QMQ - (1 + QMP) 



PM 



2x/i 



1 -1 



2^/z + i{u,Mv) 



PM 



(1 + QMP). 
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As a result, we obtain 

(57) [1-v{z-Hq)-^u]-^ 

1 - QMQ - (1 + QMP) 



2^ 



-PM 



X (1 + QMP) 



Q 



2^/z + i{u,Mv)' 
I ^ 2, 



-1 



-P + 



2^/z + i{u,Mv) 2^/z + i{u,Mv)' 



In order to obtain the explicit expression of the inverse of the operator in 
the right-hand side, we expand the operator in powers of -^/z. 
To begin with, we note that 



with 



Ko = {u,Mov)-^ {u,vR^^Kv)Po - PvR^^^u , 
where Pq = {u, Mov)-^PMo. We write Qo = 1 - -Po- This yields 

KoMoQo + QovR^o^uQo + QqMoKo] +0{e^/^+^') 



QMQ = QoMoQo- 



Therefore we have 



4i 



(58) 



where 



1-QMQ-{1 + QMP)- 



2v^ 



2^/z + ^{u, Mv) 
1 - QoMoQo + VzKi + Di 



PM 



+ 



1 

4i 



KoMoQo + QoMqKo + QovPf^uQQ 

^PoMo + QoMoPqMo 



{u,Mqv) 

and the operator Di satisfies 

\\Dih = 0{e' 

The following lemma is due to Bolle, Gesztesy and Klaus [7]: 



Lemma 17. If the operator QqMqQq has the eigenvalue +1, then the eigen- 
value must he simple and H correspondingly shows a zero-energy resonance. 

Proof. Let ipi be the eigenvector of QqMoQq with the eigenvalue +1, and 
let 

= (u, Movy \u, Mq(Pi) - A^^u(fi. 
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Prom QoMoQoipi = ipi, we have Poifi = 0, and 

¥^1 = (1 - Po)Moipi = Moipi - {u, Mov)-^v{u, ■ ■ ■ )MoVi 
= Mo^pi - {u,Mov)-\u,M^(pi)v. 

Combining this with the definition of ^, we obtain 

vip = {u,Mqv)''^{u,Mq(Pi)v — vA~^uipi 
= {u,Mov)~'^{u,M^(pi)v - Mo(pi = -ifi. 

Therefore, we get the integral equation for ip, 

if) = {u,MQv)~'^{u,MQipi) + /\~'^Vil). 

Clearly, A'0 = Vip^ i.e.. Hip = 0. In the same way as in Lemma 16, the 
uniqueness of tpi can be proved with the help of the relation (59) below. □ 

Let ipi be the eigenvector of QqMqQo with the eigenvalue +1. Then the 
left eigenvector is given by (pi = 7^o{sgnV)MQipi, where 'Nq is a normal- 
ization constant. Since one has {uA~^uipi,ipi) < 0, we can choose Kq to 
satisfy ((^i, </?i) = 1. 

In order to show that ipi is the left eigenvector, we note the following: 
Prom QqMoQocpi = cpi, we have 

= Po^i = iu,Mov)-^PMoipi 

= {u,Movy'P{Po + Qo)Mo^i 
= {u,Mov)-^pQoMoipi 
= {u,Mov)-^Pipi, 

where we have used PPq = 0. This yields 

(59) {u,ipi) = 
and 

(60) {u,Moipi) = 0. 
The latter relation yields 

(61) Po*^i = {u,Mov)-'M*P*isgnV)MoPi 

= {u, Movy'^uA^'^v ■ u{v, ■ ■ ■ ){sgnV)Mo(pi 
= {u,Movy'^uA~^v ■ u{u,Mo(pi) = 0. 
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Note that 



Q*o{sgnV)MoP 

= (sgnF)MoP - Po*(sgnF)MoP 

= (sgnF)MoP - (n, Movy^uA^'^v ■ u{v, ■ ■ ■ ){sgnV)Moviu, ■■■) 
= {sgnV)MoP - uA-'^v ■ u{u, • • • ) = 0. 

Using these relations, we obtain 

iQoMoQoT0i = QouA-\{sgaV)Moipi 
= QS(sgny)MoVi 

= Q*o{sgaV)[l - {u,Mov)-''MoP]M^ipi 
= QS(sgnF)Mo[l - {u,Mov)-^PMo]Mo^i 
= QS(sgnF)MoQoMo(^i 
= Qo(sgnF)Mo(^i = 01, 

except for the normalization of (pi. 

We define Pi = ipi [ipi ,■■■), and Qi = 1 — Pi . As an approximate inverse 
for the operator of (58), we introduce 



(62) 
where 
(63) 
and 



P = (1 + K2){1 - QoMoQoQi)-' + 



H -1 



1 



iu,Mov)-\u,Mlipi)\'^ + -\{xu, ipi)\'^ 



^{PiMoKo + PivR^o^uQo) - {u, Mov)-^PiMoPoMo 



Prom (60), we have 
This yields 

Further, from (61), we have 



{01, v) = 0. 
PiP = 0. 
PiPo = 0. 
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Using these relations, we obtain PiKo = 0. Prom these observations, we 
have 



~ 2i 
= -J^ 



liPiMoKoPi + PivR^^\Pi) - {u,Mov)~^PiMoPoMoPi 



K2Q1. 



Note that 



PiMqKqPi = -{u,Mov)-^PiMQPvRf\Pi 

= -{u,Mov)-\^i,Mov){u,vR^Q'>uipi)Pi 

= 2{u,Movy^y^*o{M^ipuu){u,xv){xu,^i)Pi, 

where we have used (Mof,(^i) = 'No{u, Mgipi), {u,v) = and {u,ipi) = 0. 
Similarly, 

PivR'^^Pi 

= {ipi,vR^Q^uipi)Pi 

= -2{ipi,xv){xu,ipi)Pi 

= ~2y^*o{Moipi,xu){xu,ipi)Pi 

= ~2J^\{xu,ipi)\'^Pi -2{u,MovY^Jil{PM'^ipi,xu){xu,ipi)Pi 
= -23^*o\{xu,<fi)\^Pi - 2{u,Mov)-^3^*o{M^ipi,u){v,xu){xu,ipi)Pi. 

Here, we have used 

(64) Mo<pi = (Qo + Po)Mo<pi = <fii + (u, Mov)-^PM^(pi. 

Moreover, wc have 

-{u,Mov)-^PiMoPoMoPi = ~{u,Mov)-\ipi,Mov){u,M^ipi)Pi 

= -KMoi;)-XlK^oVi)l'^'i- 

Putting these together, we have 



-P^Ki = Pi-K2Qi. 
2i 



Immediately, this yields 



B 



1 - QMQ - (1 + QMP) 



2a/z + i{u, Mv) 



B 



(1 - QoMoQoQi)Qi + V^Ki + Di 



PM 



I + D2, 
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where 
D2 



PiDi + + K2){1 - QoMoQoQi)-'Ki + 0{e 



1/2+5'-, 



Therefore, the desired expression of the inverse of the operator is 



1-QMQ-{1 + QMP)- 



-\ -1 



2^/z + i{u,Mv) 
Substituting this into (57), we obtain 



PM 



I + D2 



-B. 



[1 - v{z - Ho) 



1 



1 + 



-B{l + QMP) 



Q 



-P + 



P 



2^/z + i{u,Mv) 2^/z + i{u,Mv) 



Using this expression, we have 



(65) Tr {v{z - Ho)-^u[l - v{z - Hq) 



= Trv{z- Ho)-'^u 



1 



I + D2 



BQ 



+ 



2^/z + ^{u,Mv) 

2^ 
2yJz + i{u,Mv) 



Tr v{z - Ho)-'^u 



Tr v{z - Ho)-'^u 



1 



I + D2 



1 



I + D2 



B{P + QMP) 



B{P + QMP), 



where we have used PP = P. 

Consider first the first term in the right-hand side of (65). It is written 



(66) Tr v{z - Ho)~^u 



1 



I + D2 



BQ 



Tr v{z - Hq)-^u 



1 



I + D2 



{1 + K2){1 - QoMoQoQi)-^Q 



+ -^Tr v{z - Ho)-\-^PiQ 
2%^z \-\-U2 

by using the expression (62) of B. The expansion corresponding to (54) is 
given by 



v{z - Ho^u = -'-z-'/^P - iz-'/^vR^^K + 0(£-^/2+'5') 
4 8 



In the same way as in the preceding case, it is enough to treat 

Tr v{z - Hor^uQ 

for the first term in the right-hand side of (66). Using Q = 1 — P and the 
above expansion, one can show that the corresponding contribution to the 
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excess charge is vanishing. Similarly, the second term in the right-hand side 
of (66) is computed as 



(67) 



Tr v{z — Hq) 



I + D2 



PiQ 



1 



\ 8 

= -^(^i,^;it;('Wi) + o(£-'+'') 

= ^((^i,x^;)(xu,¥^) +0(£-^+'^') 



PiQ + 0(£-^+''') 



8z 



xti, (/?i)|^ + (u, Mof) "^(Mq «)(?;, xm)(xm,(/?i) J + 0(e' 



where we have used QP = 0, PiQo = Pi, {(fijv) = {u, ipi) = and (64). 

In the same way, the second term in the right-hand side of (65) is com- 
puted as 



(68) 



2^/z + i{u, Mv) 



Tr v{z - HoY^u 



1 



I + D2 



B{P + QMP) 



{u,Mqv)-^^Tt vR'^\Pi{l + QoMo)P + 0{e-^+^') 

= —^{u,Mov)-\u,xv){xu,ipi){Miipi,u) + 0{e-^+^'), 

where we have used P^ = 0, PiP = 0, PiQo = Pi, (u, v) = and {u, ipi) = 0. 
Similarly, the third term in the right-hand side of (65) is computed as 



(69) 



-Tr v{z - Ho)-'^u- 



1 



-B{P + QMP) 



= 7= \ P^-F^B{^ + QM) + 0(£-y^) 

2z2^ + i{u,Mv) I + D2 ^ ' ^ ' 



1 



2z2^^i{u, Mv) 
i 1 
2i2v^ + z(n, Mv) 
1 



Tr P 
Tr P 



1-Do 



1 



'I + D2 



2z 



3(1 + QM) + 0{e-y^) 
(l + QM) + 0(£-i+^') 



Tr PMo Tr PKiPiMq 

2i 



+ 0(e-i+^') 
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where we have used PP = P, PPi = PiP = 0, {u,v) = 0, PiQo = Pi, 
PQo = P and 

P{1 - QoMoQoQir' = P. 
The second term in the last line of (69) is calculated as 

Tr PKiPiMo 



Tr P 



1 



■^.{MoKo + vR^^^u) + 2i{u, Movy^MoPoMo 



Pi Mo 



1 1 

= - — {u,Mov)-^Tv PMoPvPf^PiMo + ^Tr PvR^qKPiMq 
+ 2i{u, Mo?;)~^Tr PMqPqMqPiMq 

where we have used PKq = 0, PiQo = -Pi, Qo-Pi = -Pi, PqP = P and 
PPq = 0. The first and the second terms in the second equality cancel each 
other out. The third term is written 

2i{u,Mov)-^Ti PMqPqMqPiMo = 2fTr PoPoAIoPiMq 

= 2iTr Po Mo Pi Mo 
= 2i{u,Mov)-\u,M^^i){ipi,Mov) 
= 2i{u,Mov)-^J<*o\{u,M^^i)f. 

Therefore, the corresponding contribution in (69) becomes 

1 



2z 



JOrQ{u,Mov)-'\iu,M^ipi)\ 



In consequence, the nonvanishing contribution from the third term in the 
right-hand side of (65) to the excess charge is 

1 T\n\r* 

-2^ + ^(^,Mo.)-2|(u,MoVi)|^. 

Combining this, (63), (67) and (68), we have the desired result, Z° = 0, for 
the excess charge at the zero energy. 



B Proof of Theorem 8 

Since a > 6, one can find a small S > such that a > 6 + S. As in 

Section 4.1, wc take z = ee^'^ with 5 < (j) <2'k — 5 with 5 > 0, and R = 
with 3/(6 + 6) < ^ < 1/2. We write Gq{x;z) for the integral kernel of 
the resolvent Rq{z) = {z — Hq)^^. We decompose Ro{z) into two parts as 
Ro{z) = Rq {z) + Rq {z) , where Rq{z) has the integral kernel Go{x; z)x^{x), 
and Rq{z) has Go{x; z)xr{x)- The integral representation of Go{x; z) yields 

(70) \Gq{x-z)x^{x)\ < Const.e'^-^/^ 

We write v = \V\^/'^ and u = \V\^/'^sg-nV . 
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Lemma 18. The following bound is valid: 

\\vR^{z)u\\^ < Const.e^+'^' 
with some positive 5' for a small s. 
Proof. Using the bound (70) , one has 

\V{x)\'/^Goix y; z)x^{x - y)\Viy)\'/hgnV{y)\^ 

< Const.|y(x)||F(y)|x^(x - y)e^''-^ 

< Const.\Vix)\{x)''+'\V{y)\{y)^+\xy^-'xU^ - y){y)-^-'e^'^-' 



4+<5 



<Const.|F(x)|(x)^+^|y(y)|(y)^+'5 
< Co^si.\V{x)\{x)^+'^\V{y)\{y)^^h^i'+^i'-^. 
Combining this with the assumption on the decay of the potential V yields 
||i;i?5'(2)tx||2 < Const.£3'^+'^'^/2-i/2_ 

Since /x > 3/(6 + (5), the desired bound holds. □ 

Next, consider vRq{z)u. To begin with, we recall the expansion of 
Gq{x\z) with respect to y/z\x\ as 



Gq{x\z) 



k{z) + — log \x\ 



1 - +^QZ^|x|^'^ 



^ oo 



1=2 



with the coefficients, q and c^, where 

1 



k{z) := 



27r 



ni 1^1, 
-y + 7 -log 2 + - log z 



with Euler's constant 7. For i >2, one has 
v{x)z^\x - yl'^^Xnix - y)u{y) 



< \v{x)z\x - yfu{y)\ 



X \z 



i-1 



\x-yf^'-'hU^-y) 



< e 

< 4e 



\V{x)\y'{x)'{x)-'\x - y\'{y)-'\V{y)\'/'{yf\ {e'-'^y-' 
\V{xr^{xf\Viy)\'/^{yf (e^-^'^)^-\ 
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Therefore the corresponding Hilbert-Schmidt norm becomes 0{£^'^^') with 
some positive 6'. This yields 



v{x)Go{x - y; z)xr{x - y)u{y) 



v{x) 



k{z) + — log \x - y\ 



1- ^k-yl^J u{y)xRix-y) 
.1+6' ^ 



+ ^zv{x)\x - y^xU^ - vHy) + )■ 

In the same way, one can replace the cutoff Xr with 1 in this right-hand 
side. Combining the result with the above Lemma 18, we obtain 

vRo{z)u = k{z)v{u, • • • ) + vA~^u — -zk{z)vR^^ u 



+ ^zvR^o^u - ^zvR^aK + 0{e 



1+5' 



where the operators, R^^ and R^\ have the integral kernels, |xp and 
|a;plog|x|, respectively; A is the Laplacian. 

Consider first the case with J (fxV{x) = {u,v) 7^ 0. We define 

P := {u, v)~^v{u, • • • ) 

and Q -.= 1- P. Write 

vRo{z)u = {u, v)k{z)P + M 

with 

(71) M:=Mo- ^Kiz)vR!^\ + ^vR^u - ^vR^u + 0{e^+^'), 

4 OTT OTT 



where we have written Mq = vA ^u. Note that 



1 - {u,v)k{z)P 



1 + 



{u,v) 



K ^{z) — (n, v) 
This yields 
l-v{z- Ho)-^u =l-{u, v)k{z)P - M 



P = Q 



K-\Z) 



{u, v) — K ^{z) 



P. 



1 - {u,v)k{z)P 



1-QM + 



{u, v) — K ^{z) 



PM 



1-{u,v)k{z)P (l-QMP) 
1 - QMQ - (1 + QMP) 



{u, v) — K ^{z) 



PM 
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Therefore the inverse is 
(72) 



[1 ~ v{z - Ho)-\] ' 



1 - QMQ - (1 + QMP)- ^''\. . PM 

{U,V) — K~'-{Z) 



X (1 + QMP) 



Q- 



{u, v) — K ^{z 

Using the expansion (71) for M, we have 



-P 



(73) 1-QMQ-{1 + QMP) 



n-\z) 



{u, v) — K ^{z) 



PM = l-QMoQ + Ki + Di, 



where 



Ki = -K{z)QvRfuQ - -^QvR^uQ + -^QvPfuQ 

'^^TrT(l + <3MoP) [PMo - jK{z)PvR^i\ 



{u, v) — K ^{z) 

+ ^{u,v)-^QvR^\PMo, 

and the operator Di satisfies ||Di II2 = 0(e| log e|~^). 

In order to obtain the inverse of the operator (73), we must consider the 
eigenvalue problem for the operator QMqQ. 

Let (f be an eigenvector of QMqQ with the eigenvalue 1, i.e., QMoQip = 
if. Clearly, one has P^p = and (n, ip) = 0. We note that 

QMoQ<f = (1 - P)Moip = Moif - PMo^ 

= Mqip — [u, v)~^{u, Mq(p)v = ip. 

The resonance state is given by [6] 

■\\) = (ti, t;)^"^(n, M{^p) — /\.~^uip. 

Actually, one can check this as follows: Using this and the above relation, 
one has 

= (n, v)^^{u, Moip)v — Moifi = —<p. 
Substituting this into the expression of ip yields 

ip = {u, v)~^ {u, Mocp) + A-^ViJj. 

Immediately, A'^ = —Vip. Besides, the asymptotics becomes 



ijj '-^ {u,v) ^{u,Mq(p) 



2tt\x\ 



J (fy yu{y)(p{y) 
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for a large |x|. Here, the leading term from the integral is vanishing because 
{u,ip) = 0. Thus, there appear at most three resonance states [6]. 

In the following, we will treat only the case that H shows three resonance 
states at the zero energy and N bound states at the zero energy because the 
other cases can be handled in a much easier way. 

Let (po,i, i = 1,2,3, be the three zero-energy resonance states satisfying 
the following conditions: 

(m,Mov9o,i) /O, 

{U, Mo(^o,2) = {U, Mo99o,3) = 0, 

and 

{XU, (/9o,2) / 0, {XU, V9o,3) / 0. 

Here x is the operator of the position. Therefore x is a two-component 
vector, and the two vectors, {xu, (^0,2) and {xu, (^0,3) are linearly independent 
so that three zero-energy resonance states appear. Let j = 1,2, ... , N, 
be the zero-energy bound states satisfying 

{u,MQ(f)j) = 0, and {xu,(f)j) = 

for j = 1,2,. . . ,N. Further we define the projection operators as follows: 

Po,i ■■= Vo,i{^o,i, ■■■) for i = 1, 2, 3 

and 

N 

Here (^o,« and are the corresponding left eigenvector for the operator 
QMoQ. They are given by 0o,i = ^o,i(sgnV)(^o,j and = Moj(sgnF)^j. 
Here, we have chosen the normalizations, Ko,i and Mqj, to satisfy 

(<^0,i, VO,k) = Sik, 

{4>j,(pe) = Sjc, 

and 

We also write Pq = Po,i + Po,2 + Po,3, P = Pq + Pi and Q = 1 - P. 

Since the two vector, (x-u, (^0,2) and (xu, (/^q^s) are linearly independent, 
there exists a matrix such that 

XI ^iji^o,j,xv) ■ {xu, ipo,k) = Si,k for i,k = 2, 3. 

j=2,3 
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Similarly, since the operator uR^^u is positive [6], there is a matrix 4>jj such 
that 

N 

XI ^i,jih^'"^o'''^M = Si,k for i, A; = 1, 2, ... , N. 

As an approximate inverse for the operator (73), we introduce 
B = {1- UKi){l - QMoQQ)-^ + U 



with 



zk(z) ^ 



J' 



Stt 



N 



+ — X] (l^k^k,ei4>e, ■■■)■ 



k,e=i 

Using the properties of the eigenvectors of QMqQ with the eigenvalue 1, one 
can show 

n^sTi = 0(1). 

This yields 

B [(1 - QMoQQ)Q + Ki + Di 

= (1 - UKi)Q + (1 - ni^i)(l - QMoQQ)-^Ki + UKi + UDi + o{e) 
= Q + UKiP + (1 - UKi){l - QMoQQ)-^Ki + UDi + 0(e). 

Similarly, the second term in the second line is calculated as 

n^iP = P + E + PD2, 

where 



N 



2= X 'Po,i^ij{^o,j,' 



)XV 



{xu, ■ ■ ■ )Po,l + X ^k^kMi, ■ ■ ■ )vRf''uPo 



k,e=i 



+ 



2tt 

{u,v) 



X ^k^kiiL ■ ■ ■ )vR^o^uPMoPo,i 



kl=l 



and the operator D2 satisfies ||i^2||2 = o{l). Therefore we have 



B 



{l-QMoQQ)Q + Ki + Di = l + E + {l-QMoQQ)~'Ki+PD3 + o{e), 



where the operator satisfies ||£'3||2 = o{l). In consequence, we obtain 



(74) 



(1 - QMoQQ)Q + K1+D1 
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1 



1 + P>4 



-TB, 



where 



L>4 = (1 - QMoQQ)-^Ki + PD'^ + o{e) 
with satisfies II-D3II2 = and 



r = i-s + s2. 



Here we have used 

(1 + H)(1-H + h2) = 1 

which is derived from = 0. 

Consequently, from (72), (73) and (74), we have 

[1 - v{z - Ho)-'u] = ^^^TB{1 + QMP) 



Q 



{u, v) — K ^{z) 



Prom this expression, we have 



(75) Tr{u(z - Ho)-^u[l - v{z - //q)" } 



= Tr |v(z - Ho)~^u 
K-\z) 



1 



l+£>4 



VBQ 



- 7 ^^^YT^T, v{z - Ho)-'u-——TB{l + QM)P . 

{U,V) - K '-{Z) { 1 + D4 J 

First, we wih show that the second term in the right-hand side of (75) 
does not contribute to the excess charge at the zero energy. To begin with, 
we note that one can obtain 



(76) 
v{z — Ho)~^u 



[u,v) 
Anz 



k{z) 1 
^ 167r 



?(3)„ 



-S'^ 



in the same way as (71). Using this expansion, one has 

- 7 T — ^'^^-1. ^ Tr v(z - Ho)-^u — i — rB(l + QM)P 

{u,v) „ ^ 1 



Anz {u, v) — K ^{z) 



Tr P- 



K-\Z) 



{u, v) — K ^{z) 



k{z) 



1 



TB{l + QM) 



167r 



Tr vRf\—^VB{l + QM)P + o(£-^). 
1 + -D4 



1 



Using the properties of the vector one can show that the contribution 
from this second term is vanishing. The first term in the right-hand side is 
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computed as 

(77) 



1 



1 + Di 



rB{l + QM) ^ = TV <^ P 



1-Da- 



Tr|prB(l + QM)| 

1 



TB{l + QM) 



l + £>4 



TB{l + QM)P 



Using Pifo^i = P<pj = 0, the first term in the right-hand side of the second 
equahty is calculated as 

TrPrS(l + QM) = TiPB{l + QM) 

= TVP(1 - QMoQQ)"^(l + QM) 
= TrP(l + QM) = TrP = 1. 

Here we have used P(l — QMqQQ)~^ = P. As to the second term, we note 
that, in the same way, we have 



and 



PD4 = PKi + o{e) = 0(| logep^), 

B{1 + QM)P = 5 [(1 + QMo) + Q{M - Mq)] P 

= (1 - UKi){l - QMoQQ)-\l + QMo)P 
+ n(l + QMo)P 



+ 



{l-UKi){l-QMoQQ)-^ + U Q{M-Mo)P 



= YIMqP + Yl{M - Mq)P + 0{l) = 0(1 log e|), 

where we have used {(pQi^M^v) = for i = 2,3 and {(f)k,MQv) = for 
k = 1,2, . . . , N . Therefore the second term in the right-hand side of the 
second equality of (77) becomes 0(1). Thus, the second term in the right- 
hand side of (75) does not contribute to the excess charge at the zero energy. 
The first term in the right-hand side of (75) is written 



(78) 



Tr <^ v{z - Hq)-^u 



1 



l + £>4 



TBQ 



1 



+ Tr\v{z- Ho) 



u 



1 



r(l - UKi){l - QMoQQ)-^Q 

ruQ 
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by using the expression of B. 

Consider the first term in the right-hand side of (78). If one can substi- 
tute the expansion (76) into the first term, it is vanishing. But this procedure 
is not justified as mentioned before. We have to treat Tr v{z — Hq)~'^uQ. It 
is calculated as 



TV v{z - Hq)-'uQ = Tr v{z - HoY^u - Tr v{z - Ho)~^uP 

{u,v) {u,v), 



Attz Attz 
= 0(|log£|), 



■Tr P + 0(|log£|) 



where we have used (76). Thus the contribution is vanishing. 

Similarly, the second term in the right-hand side of (78) is computed as 



(79) 



Tr v{z - HqY'^u 

k{z) 1 
~ 167r 



I + -D4 



rng 



Tr vRf^u 



1 



1 



+ g^Tr ^;43)^__rng + o{e-'+'' ). 

Consider first the first term in the right-hand side of (79). Clearly, the 
nonvanishing contributions come from the second and third terms in 11 of 
B. Therefore we treat only the following two: 



l-h<5'^ 



k{z) l_ 



zk{z) 



and 



k{z) 1_ 

~T~ ~ 167r 



Tr vE!§''u 



1 



N 



Q 



k,e=i 



Q. 



The first one is written 
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Thus the contribution to the excess charge becomes 2. This is nothing but 
the contribution from the two resonance states. 
The other one is written 

N 

hii^e, ■■■) 



Stt 
z 



'k{z) 



Stt 
z 



1 1 f 1 

V Tr <^ vR^^'' u rc!)fe$A 

16^J,^i I ° 1 + ^4 



k{z) 1 
^ 167r 



k,e=i 



by using the properties of the vectors 0j . The matrix elements are calculated 



as 



= {(i>i-,vRf^u 



-D4 + Di 



1 



1 + 

= -{^,, vRfu{l - QMoQQr'K^cPk) + o{\ loge\-^) 

This implies that the corresponding contribution is vanishing. 

Next, consider the second term in the right-hand side of (79). Clearly, 
the nonvanishing contribution is written 

1 ^ 1 

- ^ TV vRl^^u—^r(Pk^kA^e, ■■■)Q 



kl=l 



1 ^ N 

- ^kA^i,vRfucj)k) + o{e-^) = - +o(£-i). 



kl=l 



This implies that the corresponding contribution to the excess charge is 
equal to the number A'' of the bound states at the zero energy. 
Next, consider the case with 



(fxV{x) = {u,v) = 0. 



We write 



P = v{u,---). 

Then, we have 

l-v{z- Ho)-^u = 1 - k{z)P - M 

= [1 - k{z)P] [1 - (1 + k{z)P)M 

= [1 - k{z)P] \l-M- k{z)PM 
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Since 



{v,Mu) = {u,vA~^uv) + 0(£l0g|£|), 



the matrix element (u, Mv) is nonvanishing for a small e. Relying on this 
fact, we introduce 

P := (u^Mvy ^PM, 



and 



Q--1-P. 



These satisfy P'^ = P and QP = PQ = 0. Note that 



[1-{u,Mv)k{z)P]-^ = 1- 



{u,Mv) 



-P = Q 



{u,Mv) - K-'^{z) {u,Mv) - K-'^{z) 

Using this relation, one has 

l-v{z- Hq)-^u = [1 - k{z)P] [l - {u, Mv)k{z)P - M 

= [1 - k{z)P][1 - {u,Mv)k{z)P] 



P. 



K-\Z) 



-PM 



{u,Mv) - K-^{z) 
[1 - k{z)P][1 - (n, Mv)k{z)P]{1 - QMP) 

1 - QMQ + (1 + QMP)——^P^^—-PM 

{u,Mv) — K '■{z) 



Therefore, the inverse is given by 
(80) [1 - v{z - Ho)-\]-^ 

1 - QMQ + (1 + QMP) 



{u,Mv) - K-'^{z) 



PM 



X (1 + QMP) 



Q- 



K-\Z) 



{u,Mv) - K-^{z) 

,.-1 



[1 + k{z)P] 
-I -1 



1 - QMQ + (1 + QMP) ^ PM 

{u,Mv) — k '■{z) 



X (1 + QMP) 



Q- 



{u,Mv) - K-^{z) {u,Mv) - K-^{z) 



where we have used PP = P and QP = 0. 
We write 

M = Mo + Mi + O{e^+^') 

with 

4 ° Stt ° Stt ° 
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We have 

with pQ = {u, Mov)-^PMo and 

Kq = (u,Mov)-^ PMi- {u,Miv)Pq . 
Prom these expansions, 

QMQ = QoMoQo - KoMqQo + QoMiQo - QqMoKo + 0(£^+^'), 
where Qq = 1 — Pq. In the same way, 

K-^{Z) 



(1 + QMP) 



{u,Mv) - K-i(z) 



PM 



= (l + QoMqPo 



K [Z 



+ 
+ 



{u,Mv) - K-i(z) 
-KoMoPo + QoMiPo + QoMqKo 



PoMo 



{u, Mqv) 



PqMo 



PoMi + KqMo 



+ 0{e\loge\-^). 



Prom these observations, we have 



(81) 1-QMQ+{1+QMP) 



{u,Mv) - K,~^{z) 



PM = 1-QoMoQo+Ki+Di, 



where 



Ki = KqMqQq - QoMiQo + QoMoKo 



(l + QoMoPo 



i^-Hz) 



{u,Mv) - K-'^{z) 



PoMo 



-KoMqPq + QoMiPo + QoMoKo 

K-HZ) 



K~\z) 



(l + QoMoPo) 



{u, Mov) 



{u, Mqv) 
PoMi + KoMo 



PqMq 



and the operator Di satisfies ||Di ||2 = 0(e| log e|^^). 

In order to obtain the inverse of the operator (81), we must consider the 
eigenvalue problem for QMqQ. 

Let ip be an eigenvector of the operator with the eigenvalue 1, i.e., 
QMoQ(p = If. In the same way as in the case in one dimension, the obvious 
relation, pQip = 0, yields {u, (p) = and {u, Moip) = 0. The corresponding 
left eigenvector is given by = (sgny)Mo<^. 
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The resonance state is given by 

ijj = (n, Mov)~^{u, M^if) - A-'^u^p. 

Using the relation, 

ip = {l-Po)Mo(p = Mq^-{u,Mqv)-^PMI^p = Mo(p-{u,Mov)-\u,M^(p)v, 
we have 

vtjj = (u, Mqv)"^ (u, Mq(p)v — Mq(p = —(f. 

Therefore, 

ip = {u, Mov)~^{u, M^if) + A-Wijj. 
Immediately, A^; = Vip. Besides, the asymptotics is given by 

ip ~ {u,Mov)~^{u,M^(p) - -^^^ ■ J d^y yuiy)(p{y) 

for a large \x\. 

Similarly to the preceding case, we consider only the case that there are 
three resonance states and N bound states at the zero energy. 

Let ipo^i,i = 1,2,3, be the three zero-energy resonance states satisfying 
the following conditions: 

(n,MoVo,i) /O, 
(u, MoVo,2) = (u, MoVo,3) = 0, 

and 

{XU, (/3o,2) / 0, {XU, 990,3) / 0. 

Let (f)j be the zero-energy bound states satisfying 

{u, MqcPj) = 0, and {xu, (pj) = 
for j = 1,2,. . . ,N. Moreover, we introduce the projection operators, 

Po,i = ^o,i{'Po,i^'") fori = 1,2,3, 

and 

N 

3=1 

where (^o,i ^ind (pj are the left eigenvectors. By relying on the fact that 
{(p,(p) = {uA~^uip,(p) < 0, we have chosen the normalizations as 

{'Po,i,'Po,j) = Sij, 
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and 

{'Po,i,<Pj) = {^j,^o,i) = 0. 

We also write Pq = Po,i + Po,2 + ^0,3, P = Po + -Pi and Q = 1 - P. 
Note that 

Mo^o,i = (Po + Qo)Mo(po.i 

= {u, Mov)''^v{u, MQipo^i) + ipo^i = ipo^i 

for i = 2,3. Similarly, we have Mo0j = for j = 1,2,...,N. Relying on 
these properties, we can find two matrices, and to satisfy 

XI ^iji'Poj^xv) ■ {xu,ipo^k) = for i,k = 2,3, 

J=2,3 

and 

TV 

^^^ij{(j)j,vB!'QK(t)k) = 5ik for i,k = 1,2,. .. ,N, 
j=i 

respectively. 

As an approximate inverse for the operator (81), we introduce 



where 



p = (1 - nxi + e)(i - QoMoQqQ)-^ + n. 



2 



and 



Stt ^ 
k,e=i 



27r ^ 

= ^^^^ 1] (f)k^ui4>e,x^v){u,M^ ■■■). 

Then, we have 

P[(l - QoMoQoQ)Q + Ki + D^] 

= (1 - n^i + @)Q + (1 - UKi + e)(l - QMoQoQ)-^Ki 
+ n^i + UDi + P X 0{£) + o(e) 

= Q + ni^iP + eg + (1 - nKi)(i - QoMoQoQT^Ki + p x o(i) + o{e). 

The second term in the last line is calculated as 

UKiP = P + E + Px o(l). 
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where 

S= ^ ipo,i^ijiipoj,xv) ■ {xu,ipo,i){ipo,ir--) 

ij=2,3 

N 

k,e=i 

Using the properties of the eigenvectors of QqMqQo with the eigenvalue 1, 
we have 



(82) ni^i = ^ X] h^keik ■ ■ ■ )^vR^^KQo + 0(1) 

k,e=i 

N 

= 2TrK{z) <Pk^ke{4>e, x\){u, ■■■) + 0(1). 
k,e=i 

Therefore, 

UKi{l-QoMoQoQ)-^Ki 

N 

= 2TrK{z) Y (t>k^kl{^l,x^v){u, - ■ ■ ) 



k,e=i 



X (1 - QoMoQoQrH^ + QoMoPo)j^^^PoMo + o(l) 

{u,Mov) 



2vr ^ 



Y (Pk^ke{^e,x'^v){u, ■ ■ ■ )MqPoMq + o(l) 



27r ^ 

= I . yZ (l>k'^kii4>e,x''v){u,M^ ■ ■ - XPo,! + Q) + o{l) 

= ePo,i + eQ + o(i). 

where we have used 

{u,---){l-QoMoQoQ)-^ = {u,---) 
for getting the second equahty. Prom these observations, we have 

B[{1 - QoMoQoQ)Q + + Di] 

= 1 + E- GPo,i + (1 - QoMoQoQy^Ki +Px o(l) + o(e). 
The desired inverse is written 

[(1 - QoMoQoQ)Q + Ki + Di]-^ = 3^^r5, 
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where 
and 



D2 = {1- QoMoQoQy'Ki + P X o(l) + o{e), 



r = i-H + ePo,i + (s-ePo,i)'- 

From this, (80) and (81), we obtain 

1 



[1 - v{z - Ho)-'] 



Therefore, we have 



ii-i 



I + D2 

Q 



rB{l + QMP) 

p 



{u, Mv) - K,-^ {z) {u, Mv) - K,-^ (z) ' 



(83) TV {viz - Fo)-2n[l - viz - Ho)-']-^} 



= Triviz- Ho)-'u- 



+ D2 



TBQ 



iu,Mv) - K-^{z) 
iu,Mv) - K-^{z) 



Ti\v{z- Ho)~^u 



1 



1 + D; 



-TB{l + QM)P 



Tr <! viz - iJo)-\^-^rP(l + QM)P\ , 



where we have used PP = P. 

Consider first the first term in the right-hand side. Substituting the 
expression of B into it, one has 



(84) Ihilviz- Ho)~^u 



I + D2 



TBQ 



TT\viz-Ho)-^u 



1 



I + D2 



r(l - UKi + G)(l - QoMoQoQ)-^Q 



+ Tr{viz- HqY^u 
Note that 



1 



viz- Hq)-^u = -— + 



I + D2 

niz) 



TUQ 



1 

16^ 



Since QP = 0, one formally obtains that the first term in the right-hand 
side of (84) does not contribute to the excess charge. Let us justify this fact. 
Because of the same reason as in the previous cases, it is sufficient to show 
that the contribution from Tr viz — Hq)~^uQ is vanishing. Using the above 
expansion and (u, v) = 0, we have 

Tr viz - Hq)-^uQ = Tr i;(z - Hq)-^u -Trviz- HqY^uP 

= 0+-^TrPP + O(|log£|) = 0(|log£|). 
47rz 
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Thus the statement holds. 

Similarly, the second term in the right-hand side of (84) is written 



(85) 



|w(z - 


- i^o)-' 




1 


4 


IGtt 


Stt 





1 



1 + D2 



TUQ 



I + D2 



1 + D. 

rngj + o(£ 



-rng 



As to the first term in the right-hand side, the nonvanishing contributions 
come from the two summations in 11 of 5. The first one is computed as 



k{z) _ 1 
^ 167r 



zk{z) 



Q 



1 2 

= E ^ij{'Po,j,vR^o^uT<fo,i) + o{e-^) = - + o{e-^). 



This leads to the excess charge 2. As to the other one, we have 

1 ^ 

TV vR^^K-——r V (Pk^kii^i, ■■■)Q 
J- + iJi /—"^ . 



1 ^ 

= IV vR^^K^^-^r 'Pk^keih, ■■■) + 0{e\ \oge\ 



k,e=i 



Y ^ke{4>e,vRo^u _^_^^^ (l)k) + 0{e\ loge\), 



k,l=l 



except for the pref actor which is 0(£ ^| log£|). The matrix elements in the 
last line is calculated as 



1 



= {4>e, vRqK 



1-D2 + Di 



1 



I + D2 



4>k) 



,vR^^\D2ci)k) + o{\\ogs\-^] 



= o(\ log el 



where we have used the properties of (pj. As a result, the contribution is 
vanishing. 
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Next, consider the second term in the right-hand side of (85). Clearly, 
the nonvanishing contribution is 



1 18^ 

—TV vR^^\j-—-T x—Yl ^k^ke{4>e, ■■■)Q 



1 ^ - . 

- Yl ^kei^i^ vli^K(t)k) + o{e) 



z 

k,i=l 

= - + ois). 

The corresponding contribution to the excess charge is equal to the number 
N of the bound states at the zero energy. 

Finally, wc will show that the second and third terms in the right-hand 
side of (83) do not contribute to the excess charge. The second term is 
written 

Tr v{z - Ho)-^u—^rB{l + QM)P 
1 + D2 



k{z) _ 1 



Tr vR^^^u—^TB{l + QM)P 
1 + D2 



+ ^Tr vR^^\-^rB{l + QM)P, 

except for the prefactor. In the same way as in the above, one can show 
that these yield a vanishing contribution. 

Similarly, the third terms in the right-hand side of (83) is written 

Tr v{z - Ho)-^u—^TB{l + QM)P 

^ -Tr P ^ VB{1 + QM)P + o{£-^) 



Attz I + D2 

except for the prefactor which is 0(| logej^^). Therefore, it is sufficient to 
show that the trace of the operator in the first term in the right-hand side 
is 0(1). The trace of the operator is computed as 

(86) TV P-^—rB{l + QM)P 
i + D2 

= TV P — - — TB + TV P — - — TBQM 
I + D2 I + D2 ^ 

= TV P — ^--TB + TV PTBQM - TV PD2 — ^-—FBQM. 

1 -|- Z?2 1 + D2 
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Note that 

BP = [(1 - UKi + e)(l - QoMoQoQ)-^ + u\P 
= (1 - UKi + e)(l - QoMoQoQ)-^P 

= (i-nKi + e)p = o(i), 

where we have used (82). This imphes that the first term in the right-hand 
side of the second equality of (86) is 0(1). The second term in the right-hand 
side of the second equality of (86) is calculated as 

TrPTBQM = TrPBQM = TrP(l - QoMoQoQ)-^QM = 0(1). 

As to the third term, we obtain 

{u,--- )D2 = {u,--- )K, + o{e) = 0{\ log ej"^), 

and 

BQMv = B[Qo - Ko][Mo + 0(e| \oge\)]v + 0{£^') 
= BQqMov - BKqMov + 0(| logel) 

= (1 - UKi + G)(l - QoMoQoQy^QoMov + UMqv + 0(| loge|) 
= 0{\\oge\). 

These imply that the third term in the right-hand side of the second equality 
of (86) is 0(1), too. 

C High Energy Asymptotics in Three Dimensions 

In order to give a proof of Theorem 10, we take the contour in the expression 
(13) of the SSF to be a circle which is parametrized by ^; = Ae*"^ with a large 
A and the angle cf). 

To begin with, we note that 

(87) ^ '— = ^v- ' 



z — H z — Hn z — H z — Ho 



1 .V^ + ^V^F- 1 



z — Hq z — Hq z — H z — Hq z — Hq 



z — Hq z — Ho z — Ho z — Ho z — Ho 



z — H z — Ho z — Ho z — Ho 
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For z ^ cr{Ho), one has 



1 1 I f 1 f c giVz\y-x\ 



(88) TV — ^—V — 3— = [ (fx [ (fy 

^ ' z-Ho z-Ho (47r)2y J ^ 



V{y)- 



Combining this with 



one obtains 



x-y\ \y-x\ 

1 f n f n (,2i^\x-y\ 

^TTZ J ^Z TTZ 

1 1 



(89) /" Tr — ^—V — ^— = / (fyV(y). 

Thus the leading term have been obtained. 
For z ^ a{Ho), 



z — Hq z — Hq Z — Hq 

1 / 1 ^2 

TrV V 



z-Hq \z-Ho^ 



.Qi-Vzly-xl 



{'i-n-y J J \x — y\ 2iyfz 

=2iv^|g| 

where we have changed the variables to p = x, q = y — x for getting the last 
equality. Using the polar coordinate, q = rco with r = \q\ and co G S^, and 
integrating by parts, one has 

/ d-^q — -- — V(p + q)= du drr^ V(p + ru) 

J \q\ Js2 Jq r 

oo 




= / dio^^e^'^'^ -rVip + rio 

- / duj—^ / dre^'^'''^rV{p + 
J§2 Jq dr 



70 



The first term in the right-hand side of the third equality becomes —TrV{p) / z 
from the assumption (11) for the potential V. Therefore the contribution 
corresponding to the total phase shift ^(A) is 



(90) 

by using 



1 



167r2AV2 



d^xV{xf 



liri J 



Let us show that the contribution of the second term becomes o(A ^/^). It 
is sufficient to show that 



(91) / dcp du dr 

Jo Js^ Jo 

for A ~ DO. Write 







52 d 52 

The left-hand side of (91) can be written 

f27r 



(92) 



/ dcj) dui drexp -2AVVsin((/)/2) |/(p,g)| 
Jo Js"^ Jo '- -I 

= del) duj drexp -2A^/Vsin((/)/2) \f{p,q)\ 
Jo iS2 Jo L J 

2Ai/Vsin((^/2)l |/(p,g)| 



27r 



+ / d(p du drexp 



+ 



27r 



do; 



dr exp 



-2AVVsin(0/2)l |/(p,g)|. 



Here we take e = A^^/^ and take a to be some positive cutoff. Clearly the 
first term in the right-hand side becomes 0(A~^/^). The second term can 
be estimated as 



/ d(j) duj drexp -2A^/Vsin((/)/2) \fip,q)\ 

Jo JS2 Je L J 

dr J d^exp[— 2A^/^r sin(^] 

dr / #exp[-4A^/V^/7r] 
<47r2sup|/(p,g)|A-V2log(aAV2). 
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The third term can be written 

r2-K 



d(p / dio 


2ir 



dr exp 



-2A^/Vsin(<^/2) \f{p,q)\ 



= f^dcP [ d-\^x^(9)cxp [-2Ai/2|g|sin(0/2)l \fip,q)\ 
Jo J 191 L 



\x - y\ 



-2AV2|a;-y|sin(0/2)l \f{x,y-x)\. 



This is also vanishing as A — t- oo by using the same estimate and the as- 
sumption (11) on the potential V . Consequently, we obtain 



(93) 



— I dz Tr — ^—V — ^—V — ^— 

2m J z — Hq z — Hq z — Hq 



for a large A. Therefore it is enough to show that the contribution from the 
third term in the right-hand side of the third equality of (87) gives o(A~^/^). 
Note that 



(94) 



Tr 



1 



:V- 



1 



-V- 



1 



-V- 



1 



Z — H Z — Hq Z — Hq Z — Hq 
1 1 1 1 



= Tr 
+ Tr 
+ Tr 



-V- 

Z — Hq Z — Hq Z — Hq Z — Hq 

1 .V. 1 ■■ 1 ■• 1 



-V- 

) ■ 

-V- 



-V- 
I 

-V- 



-V- 



Z — Hq Z — Hq Z — Hq Z — Hq Z — Hi 



:V- 



-V- 



-v^v^v- 



Z — H Z — Hq Z — Hq Z — Hq Z — Hq Z — Hq 



Using u and v given by (15) and the relation (16), the third term in this 
right-hand side is written 



Tr v{z - HqY'^uW + v{z - H)-^u\v 



1 



1 



Z — Hq Z — Hq Z — Hq Z — Hi 



■u. 



The operator uRq{X it 0)v is Hilbert-Schmidt class, and the following two 
bounds hold: 

\\vRq{z)u\\ < Const.|z|"^/2 



and 



< Const. Ul ^ 



for a large \z\. (See Theorem 8.1 of [23].) Combining the bound about 
the resolvent Rq{z) with the relation (17), one notices that v{z — H)~^u is 
uniformly bounded for a large A. From these observations, the contribution 
of the third term can be estimated as 0(A~^). 
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Since both of the first and the second terms in the right-hand side of 
(94) can be estimated in the same way, we treat only the first term. Note 
that 

1 1 1 1 



Tr 



-V- 



-V- 



-V- 



z — H.Q z — H.Q z — Hq z — Ho 



Tr V 



1 



-V 



1 



-V 



Z — Hq Z — Hq \Z — Hq 

^ j (fxi J SX2 J d^X3V{xi 



J^\xi-X2\ 



-Vix2)- 



i^\x2-xz\ 



\X\—Xl\ \X2 — X-},\ 

Z 

Following [11], we change the variables as 

X\ = Xl 

X2 = Xi+ rOJ2 

X3 = X2 + roj^ = x\+ r{ui2 + wa) 

with r > and with the constraint |t(J2p + |w3p = 1. Clearly, the constraint 
implies u := {uj2,cos) G S^. Then we have 



J (fxi j Sx2 j (fx^Vixi) 

= [ (fxiVixi) [ 
J Js 



pi^/z\xi-X2\ pi\fz\x2-xo,\ 

v{x2)- -y(a;3)e^^l^3-^il 



Xl - X2 



\X2 - X3\ 



du- 

5 |W2||W3| Jo 



drr exp[i\/zr0] 



X V{xi + ruj2)V{xi + r{uj2 + ws)), 
where we have written B = |aj2| + \oj-i\ + \oj2 + '^sl- We write 

/(r) = r^V{xx + ru2)V{xi + r{uj2 + u^)), 
and evaluate the integral about r as 



^ ^ ^^/^e 



fir] 



1 



i\fzQ Jq 



d 



dre'^''®—f{ 



dr' 



+ 



2 



«9 



dr 



fir] 



1 \ ^ foo q2 



by integrating by parts repeatedly. In the same way as in the above, one 
can show that the contribution to the total phase shift becomes o(A^^/^). 
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